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ABSTRACT 


The effects of oscillating flow on the pressure force normal to the 
chord of a symmetrical airfoil were investigated experimentally employ¬ 
ing a remote pressure transducer to measure the instantaneous pressure 
distribution. 

An open circuit wind tunnel having a set of rotating shutter blades 
located down stream of the test section was used to produce the oscil¬ 
lating flow. Electrical signals analogous to the free stream velocity and 
surface pressure were recorded simultaneously on separate tracks of a 
magnetic tape. The recorded data were converted to digital representation, 
and numerical techniques utilized to evaluate the spectral composition of 
the measured pressure distribution, from which the normal force was 
calculated. 

It was found that the magnitude of the total normal force at high 
angles of attack is significantly greater in oscillating flow than in 
steady flow and is frequency dependent; while at low angles of attack 
no significant differences were observed. Moreover it was found that 
higher order harmonics of the fundamental free-stream frequency constitute 
a significant fraction of the normal force, and these fractions are also 
frequency dependent. The observed results are not adequately predicted 
by quasi-steady aerodynamic analysis. 



1 


ARY 

L POSTGRADUATE SCHOOE 
EREY, CALIF. 93940 


DUDLEY KNOX LIBRARY 
HAY'LL POSTGRADUATE SCHOOL 
MONTEREY. CA 93943-5101 


TABLE OF CONTENTS 


I. INTRODUCTION.15 

II. BASIC APPROACH.16 

A. ALTERNATIVES AVAILABLE FOR MEASUREMENT.16 

B. BRIEF DESCRIPTION OF EXPERIMENTAL METHOD . 17 

III. EXPERIMENTAL EQUIPMENT . 18 

A. WIND TUNNEL .18 

1. General Description . 18 

2. Rotating Shutter Valve . 18 

3. Test Section.22 

B. MODEL .22 

C. PRESSURE MEASURING SYSTEM .27 

1. Pressure Transducers . 27 

2. Power Supplies and Amplifier . 31 

D. HOT WIRE ANAMOMETER .31 

E. TAPE RECORDER .37 

F. MISCELLANEOUS .37 

IV. PROCEDURES .39 

V. DATA ANALYSIS .40 

A. ANALOGUE TO DIGITAL CONVERSION . 40 

B. TAPE CONVERSION .43 

C. ADDITIONAL ASSUMPTIONS . 43 

D. NON-DIMENS lONALIZATION.44 

1. Pressure Coefficients . 44 

2. Norma] Force Coefficients . 44 

3. Center of Pressure .45 


3 




























vr/ ya* •>' 

-ta fCi^- ■ ” 

i c: '■-••■ 

E. COMPUTATIONAL PROCEDURES . 45 

VI, RESULTS AND DISCUSSION .48 

A. CONDITIONS INVESTIGATED . 48 

B. STEADY FLOW RESULTS.48 

C. UNSTEADY FLOW RESULTS.49 

1. Run 5 .49 

2. Run 7 .51 

D. COMPARISON OF RESULTS AT DIFFERENT FREQUENCIES . 52 

1. Normal Force and Center of Pressure . 52 

2. Pressure Distributions . 53 

E. COMPARISON OF UNSTEADY FLOW RESULTS WITH 

STEADY FLOW ANALYSIS .54 

F. SUMMARY OF RESULTS . 55 

VII. CONCLUSIONS.57 

APPENDIX A - Transducer Calibration . 128 

APPENDIX B - Analogue to Digital Conversion . 132 

APPENDIX C - Numerical Analysis of Data.137 

APPENDIX D - Evaluation of Experimental Techniques . 140 

COMPUTER PROGRAM LISTINGS . 144 

A. DATA ANALYSIS .144 

B. TAPE CONVERSION.160 

C. DISPLAY.162 

LIST OF REFERENCES .170 

INITIAL DISTRIBUTION LIST.171 

FORM DD 1473 173 


4 



























LIST OF TABLES 


Table 


Page 

I 

Typical Tape Recorder Calibration 

37 

II 

Typical Recorded D, C. Voltage and Corresponding 
Digital Mean Value 

43 

III 

Experimental Conditions for which Data was Taken 

58 

IV 

Chordwise Variation in Amplitude Ratio, Run 5 

142 

V 

Chordwise Variation in Amplitude Ratio, Run 7 

143 


5 



I 

I 


I 




Fi; 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 


LIST OF ILLUSTRATIONS 


Page 


Plan View of Wind Tunnel 19 

Photograph of Rotating Shutter Valve 21 

Photograph of Wind Tunnel Test Section 23 

Typical Wind Tunnel Test Section Velocity Profile 24 

Overall Photographic View of Wind Tunnel 25 

Schematic Diagram of Airfoil Model Showing Pressure 26 

Tap Location 

Photograph of Airfoil Model Mounted in Tunnel 28 

Sectional Drawing of Pressure Transducer 29 

Photograph of Transducer Components 30 

Transducer Gain 32 

Transducer Phase Angle 33 

Schematic Diagram of Pressure Measuring System 34 

Typical Pressure Measuring System Static Calibration 35 
Curve 

Photograph of Pressure Measuring System Installed 36 

Schematic Diagram of Data Acquisition System 38 

Filter Gain 41 

Filter Phase Angle 42 

Abbreviated Flow Chart of Computer Program Used to 47 

Analyze Data 

Cj^ vs. a, Steady Flow 59 

Mean Pressure Distribution on Airfoil, Run 2 60 

a = 5 Degrees, Steady Flow 

Moan Pressure Distribution on Airfoil, Run 3 61 

a = 10 Degrees, Steady Flow 


7 



Fi; 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 


Page 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 


Mean Pressure Distribution on Airfoil 
a = 15 Degrees, Steady Flow 

PSD of Unsteady Pressure Coefficient, 

§ = 0.0 


PSD 
§ = 

PSD 
§ = 

PSD 

= 

PSD 
§ = 


of Unsteady Pressure Coefficient, 
0.250, Upper Surface 

of Unsteady Pressure Coefficient, 

0.0 


of Unsteady Pressure Coefficient, 
0.150, Lower Surface 

of Unsteady Pressure Coefficient, 

0.0 


, Run 4 

Run 1 

Run 1 

Run 2 

Run 2 

Run 3 


PSD of Unsteady Pressure Coefficient, Run 3 
i = 0.450, Upper Surface 

PSD of Unsteady Pressure Coefficient, Run 3 
§ = 0.5, Lower Surface 

PSD of Unsteady Pressure Coefficient, Run 4 

§ = 0.0 


PSD of Unsteady Pressure Coefficient, Run 4 
? = 0.500, Upper Surface 

PSD of Unsteady Pressure Coefficient, Run 4 
§ = 0.125, Lower Surface 


Mean Square of Unsteady Pressure Coefficient vs. 
Run 3 


Mean Square of Unsteady Pressure Coefficient vs. 
Run 4 


Typical Unfiltered PSD, Steady Flow 

Velocity and Pressure Analogue Waveforms, Run 5? 
f^ = 94 Hertz, a = 15 Degrees, g = 0.085 

Velocity and Pressure Analogue Waveforms, Run 5? 
(continued) 

PSD of Unsteady Velocity, Run 5 

PSD of Unsteady Pressure Coefficient, Run 5 
^ “ 0.0, Upper Surface 


8 



Fi 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 


Page 


PSD 
§ = 

of Unsteady Pressure 
0.025, Upper Surface 

Coefficient, 

Run 

5 

80 

PSD 
§ = 

of Unsteady Pressure 
0.075, Upper Surface 

Coefficient, 

Run 

5 

81 

PSD 
§ = 

of Unsteady Pressure 
0.125, Upper Surface 

Coe fficient, 

Run 

5 

82 

PSD 
§ = 

of Unsteady Pressure 
0.175, Upper Surface 

Coefficient, 

Run 

5 

83 

PSD 
§ = 

of Unsteady Pressure 
0.225, Upper Surface 

Coefficient, 

Run 

5 

84 

PSD 
? = 

of Unsteady Pressure 
0.300, Upper Surface 

Coefficient, 

Run 

5 

85 

PSD 

5 = 

of Unsteady Pressure 
0.350, Upper Surface 

Coefficient, 

Run 

5 

86 

PSD 
§ = 

of Unsteady Pressure 
0.400, Upper Surface 

Coefficient, 

Run 

5 

87 

PSD 
§ = 

of Unsteady Pressure 
0.450, Upper Surface 

Coefficient, 

Run 

5 

88 

PSD 
§ = 

of Unsteady Pressure 
0.500, Upper Surface 

Coefficient, 

Run 

5 

89 

PSD 
§ = 

of Unsteady Pressure 
0.600, Upper Surface 

Coefficient, 

Run 

5 

90 

PSD 
§ = 

of Unsteady Pressure 
0.700, Upper Surface 

Coefficient, 

Run 

5 

91 

PSD 
? = 

of Unsteady Pressure 
0.750, Upper Surface 

Coefficient, 

Run 

5 

92 

PSD 
? = 

of Unsteady Pressure 
0.0, Lower Surface 

Coefficient, 

Run 

5 

93 

PSD 
§ = 

of Unsteady Pressure 
0.050, Lower Surface 

Coefficient, 

Run 

5 

94 

PSD 
§ = 

of Unsteady Pressure 
0.100, Lower Surface 

Coefficient, 

Run 

5 

95 

PSD 
§ = 

of Unsteady Pressure 
0.150, Lower Surface 

Coefficient, 

Run 

5 

96 

PSD 
F = 

of Unsteady Pressure 
0.200, Lower Surface 

Coe ff icient , 

Run 

5 

97 


9 


Figure 


Page 

58 

PSD of Unsteady Pressure Coefficient, Run 5 

5 = 0.250, Lower Surface 

98 

59 

PSD of Unsteady Pressure Coefficient, Run 5 

5 = 0.350, Lower Surface 

99 

60 

PSD of Unsteady Pressure Coefficient, Run 5 
§ = 0.450, Lower Surface 

100 

61 

PSD of Unsteady Pressure Coefficient, Run 5 
§ = 0.600, Lower Surface 

101 

62 

PSD of Unsteady Pressure Coefficient, Run 5 
§ = 0.750, Lower Surface 

102 

63 

PSD of Unsteady Pressure Coefficient, Run 5 

5 = 0.850, Lower Surface 

103 

64 

PSD of Unsteady Pressure Coefficient, Run 5 
§ = 0.950, Lower Surface 

104 

65 

Mean Square of Unsteady Pressure Coefficient vs. 

Run 5 

105 

66 

Unsteady Normal Force Coefficient and Unsteady 
Velocity vs. Time, Run 5 

106 

67 

PSD of Unsteady Normal Force, Run 5 

107 

68 

Mean Presuure Distribution on Airfoil, Run 5, 
a = 15 Degrees, f^ = 94 Hertz, e = 0.085 

108 

69 

Time Dependent Pressure Distributions on Airfoil, 

Run 5 

109 

70 

Time Dependent Pressure Distributions on Airfoil, 

Run 5 (continued) 

110 

71 

PSD of Unsteady Pressure at the Center of Pressure, 

Run 5 

111 

72 

Velocity and Pressure Analogue Waveforms, Run 7, 
f^ = 11 Hertz, a =* 15 Degrees, e = 0.072 

112 

73 

Velocity and Pressure Analogue Waveforms, Run 7 
(continued) 

113 

74 

PSD of Unsteady Velocity, Run 7 

114 

75 

PSD of Unsteady Pressure Coefficient, Run 7 

5 = 0.0, Upper Surface 

115 

76 

PSD of Unsteady Pressure Coefficient, Run 7 

116 


§ =• 0.100, Upper Surface 


10 


Figure Page 

77 PSD of Unsteady Pressure Coefficient, Run 7 117 

5 = 0.500, Upper Surface 

78 PSD of Unsteady Pressure Coefficient, Run 7 118 

5 = 0.900, Upper Surface 

79 PSD of Unsteady Pressure Coefficient, Run 7 119 

5 = 0.100, Lower Surface 

80 PSD of Unsteady Pressure Coefficient, Run 7 120 

§ = 0.500, Lower Surface 

81 Mean Square of Unsteady Pressure Coefficient vs. 121 

Run 7 

82 Unsteady Velocity vs. Time, Run 7 122 

83 Unsteady Normal Force Coefficient vs. Time, Run 7 123 

84 PSD of Unsteady Normal Force, Run 7 124 

85 PSD of Unsteady Pressure at the Center of Pressure, 125 

Run 7 

86 Mean Pressure Distribution on Airfoil, Run 7, 126 

a = 15 Degrees, f^ = 11 Hertz, e = 0.072 

87 ' Mean Normal Force Coefficient vs. Frequency 127 

a = 15 Degrees, U = 100 feet/sec, e = 0.08 

88 Mean Center of I^essure vs. Frequency 127 

a = 15 Degrees, U = 100 feet/sec, e = 0.08 

89 Pressure Transducer Static Calibration Instrumentation 130 

90 Pressure Transducer Dynamic Calibration Instrumentation 131 

91 Logic and Analogue Circuitry for Steady Flow Analogue 135 

to Digital Conversion 

92 Logic and Analogue Circuitry for Unsteady Flow 136 

Analogue to Digital Conversion 


11 


1 



m 





Is. 






I 


I * 

* 

m 








m 


§ 



LIST OF SYMBOLS 


Symbol 


Definition 


N 


"N 


c 

G(f) 


H(f) 

H, 


N 

N 


P 

P 


U 

U 


c 


Normal Force Coefficient, - — 

l/apU'^c 

N 

Mean Normal Force Coefficient, -^— 

l/2pU'^c 

p 

Pressure Coefficient, -^ 

l/2pU 

p 

Mean Pressure Coefficient, - ^ 

l/2pU 

Nyquist Cutoff Frequency 

Complex Transducer Transfer Function 

Complex Value of G(f) at kth Discrete Frequency 
Increment 

Complex Filter Transfer Function 

Complex Value of H(f) at kth Discrete Frequency 
Increment 

Aerodynamic Force Normal to Airfoil Chord, Nfn 
Mean Aerodynamic Force Normal to Airfoil Chord 
Absolute Pressure, Pfp 
Mean Absolute Pressure 
Atmospheric Pressure 

Directly Measured Calibration Chamber Pressure 
Indirectly Measured Calibration Chamber Pressure 
Time 

Free Stream Velocity, IH-u 
Mean Free Stream Velocity 
Chord Length 


Units 


Hertz 


Ib/ft 

Ib/ft 

Ib/ft^ 

Ib/ft^ 

Ib/ft^ 

Ib/ft^ 

Ib/ft^ 

sec 

ft/sec 
ft/sec 
ft 


13 







Symbol Definition 

c Unsteady Normal Force Coefficient, 


l/2pU^c 


Unsteady Pressure Coefficient, 


l/2pU^c 


f 

^0 

n 

P 

q 

u 

-2 

u 

X 

y 

a 

G 

§ 

U 

P 


Frequency 

Fundamental Oscillatory Frequency of Free 
Stream Velocity 

Unsteady Normal Force 

Unsteady Pressure 

Free Stream Stagnation Pressure Defined with Mean 
Free Stream Velocity, l/2pU^ 

Unsteady Velocity 

Mean Square of Unsteady Velocity 

Chordwise Coordinate of Point on Airfoil Surface 
Measured from Leading Edge 

Coordinate of Point on Surface of Airfoil Measured 
Normal to Chord 


Angle of Attack 

Amplitude Ratio of Free Stream Oscillation, 



Nondimensional Chordwise Coordinate, x/c 
Nondimensional Coordinate Normal to Chord, y/c 
Air Density 


Subscripts 


I Refers to Lower Surface of Airfoil 

u Refers to Upper Surface of Airfoil 

Superscript 

* Denotes Complex Conjugate 


Units 


hertz 

hertz 

Ib/ft 

Ib/ft^ 

Ib/ft^ 

ft/sec 
ft /sec 
ft 

ft 

degrees 


slug/ft^ 


14 








I. INTRODUCTION 


The aerodynamic forces associated with unsteady flows have attracted 
great interest in recent times. With the rapid development of the heli¬ 
copter and advances in gas turbine technology, more detailed knowledge 
of unsteady phenomena is required than can be obtained from quasi-steady 
flow analysis. 

Experimental investigations of unsteady flow problems have been 
hampered by lack of adequate facilities, i.e. lack of wind tunnels that 
are capable of establishing a controlled unsteady flow. Consequently 
the majority of previous investigations have been concerned with the 
starting process or have utilized programmed model motion rather than 
an oscillating flow [Ref. 1, 2, and 3]. Analytical investigations of 
the problem have been restricted to linearized solutions that are not 
applicable for airfoils at large angles of attack [Ref. 4 and 5]. 

The objective of the work reported here was to measure the time- 
dependent pressure distribution on an airfoil at constant angle of 
attack in a controlled oscillating stream and from these measurements 
to calculate the time-dependent forces acting on the airfoil surface. 
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II. BASIC APPROACH 


A. ALTERNATIVES AVAIIABLE FOR MEASUREMENT 

There are two obvious methods one might use to measure the instan¬ 
taneous pressure distribution on the surface of an airfoil model in an 
unsteady flow. The most direct technique would be to mount many trans¬ 
ducers in the model in such a manner that the pressure at discrete 
points on the surface could be measured directly. The difficulties 
involved in applying this method are the cost of the miniature trans¬ 
ducers required and the possibility that the dynamics of the test 
installation would influence the response of the transducers. 

A less direct method would be to connect remote transducers to 
pressure taps on the surface of the airfoil with tubing. This would 
eliminate the dynamic interference likely to be encountered with the 
more direct method, but would add to the complexity of the system 
because of the transfer functions associated with the tubing. 

In the case of oscillating flows, however, both the free stream 
velocity and the pressure at a point on the surface of the airfoil are 
periodic functions of time having some fixed phase relation between 
them. Knowledge of the phase relations existing between the free 
stream velocity and the pressure at various points on the surface of 
the airfoil, assuming that this phase relation remains constant, would 
allow the pressure distribution to be reconstructed from temporal 
records of the pressures and velocities taken over different intervals 
of time. A remote transducer then could be utilized to measure pres¬ 
sure and, if the leads from each tap were identical, only a single 
transfer function would have to be known. 
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B. BRIEF DESCRIPTION OF EXPERIMENTAL METHOD 


The technique utilized in this investigation was based on the last 
alternative described above. Electrical signals analogous to the free 
stream velocity and surface pressure at a point were recorded simul¬ 
taneously on separate tracks of magnetic tape. The pressure analogues 
recorded were the electrical output of a single remote transducer con¬ 
nected in turn by tubing to each of the pressure taps on the surface of 
the airfoil. The velocity analogues recorded were the output of a 
linearized hot-wire anemometer. 

The transfer function of the tubing was experimentally measured and 
applied in digital form to the Fourier Transform of the digitized 
pressure analogue in order to obtain the Fourier Transform of the 
pressure at the surface of the airfoil. The inverse of the pressure 
Fourier Transform was then taken to obtain the pressure-time relation. 


17 


Ill, EXPERIMENTAL EQUIPMENT 


A. WIND TUNNEL 

1. General Description 

The experimental work was conducted in the low-speed, oscillating 
flow wind tunnel located in the Aeronautics Laboratories of the Naval 
Postgraduate School. This wind tunnel is of open circuit design, with 
a 24-inch square by 223-inch long test section. A plan view of the 
tunnel is presented in Figure 1. The tunnel inlet is eight feet square, 
resulting in a 16:1 contraction ratio. Three high solidity screens 
located in the inlet section just upstream of the nozzle produce measured 
free stream turbulence intensities of 0.261 to 0.413 percent for the 
velocities encountered in the present work. 

The wind tunnel drive consists of two Joy Axivane Fans in series, 
each of which has an internal, 100 horsepower, direct connected, 1750 rpm 
motor. The fan blades are internally adjustable through a pitch range of 
25 to 55 degrees, providing a wide operating base. Two sets of variable 
inlet vanes, located immediately upstream of each fan, are externally 
operated to provide control of test section velocity. These vanes are 
of multileaf design, and preswirl the air in the direction of fan 
rotation to reduce fan capacity. The total range of tunnel velocity is 
from 10 to 250 feet per second. 

2. Rotating Shutter Valve 

Two fundamental methods of creating an oscillating flow environ¬ 
ment have been employed in the past. Both Nickerson [Ref. 6] and Hori 
[Ref. 7] introduced oscillations by oscillating their models in a steady 
flow environment. This method severely restricts the range of attainable 
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PLAN VIEW OF WIND TUNNEL 























frequencies because of mechanical complications, and also introduces 
measurement difficulties. The other approach is to actually oscillate 
the flow over a stationary model. Hill [Ref. 8] used a sliding shutter 
to impose oscillations on the free stream but was restricted by mechanical 
limitations to low frequencies. 

The most successful method of obtaining an oscillating flow with 
large ranges of frequency and amplitude was that employed by Karlsson 
[Ref. 9], and later by Miller [Ref. 10] in his investigation of trans¬ 
ition. A rotating shutter valve, immediately downstream of the test 
section, is employed to superimpose a periodic variation of velocity on 
the mean flow. The method used in the present investigation is virtually 
identical to that employed by Miller. The shutter valve consists of four 
horizontal steel shafts equally spaced across the test section. The 
shafts are slotted to accommodate flat blades of various widths, forming 
a set of four butterfly valves spanning the test section. Figure 2 is a 
photograph of the shutter valve. Each blade is driven from its immediate 
neighbor by means of a timing belt and pulley arrangement. The bottom 
shaft is driven by a five horsepower variable speed electric motor, 
through a timing belt and pulley. An intermediate shaft between the 
motor and shutter valve permits a wide variety of pulley ratios. This 
drive arrangement provides a frequency range of two cycles per second to 
the first critical frequency of 933 cycles per second. The electric 
motor presently in use, however, restricts the oscillation frequency to 
a maximum of 240 cycles per second. The amplitude of oscillation is con¬ 
trolled by blade width. Test section closure may be varied from 25 to 100 
per cent. The resulting amplitude of oscillation of test section velocity 
is a function of frequency, mean velocity and pressure gradient. In this 
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GjRE z photograph of rotating shutter valve 


investigation, blades producing 50.0 per cent closure were used, resulting 
in an amplitude range of from'3 to 40 per cent of the local mean free 
stream velocity. 

3. Test Section 

The wind tunnel test section is shown in Figure 3. Continuous 
pieces of two-inch thick aluminum, 24 inches wide and 223 inches long, 
form the upper and lower test section walls. Each of the side walls 
consists of three two-inch thick panels of stress-relieved Lucite. For 
this investigation, the central side wall panel on the opposite side of 
the tunnel from the control console was replaced with two-inch thick 
plywood to facilitate the mounting of instrumentation. The Lucite panels 
on the console side of the test section are hinged and may be raised 
hydraulically, providing access to the test section. The heavy con¬ 
struction of the test section is intended to minimize deflections 
induced by rapid changes in static pressure. 

Figure 4 is a typical test section velocity profile. Velocity 
variation is less than one per cent from the mean to within three inches 
of any wall. 

Figure 5 is an overall photographic view of the wind tunnel. 

B. MODEL 

The airfoil model used in this investigation was that employed by 
Allen [Ref. 11]. It is a NACA 63-010 section modified by straight line 
fairing from 60 per cent chord to the trailing edge. The model has a 
24 inch span and a constant 6 inch chord. Twenty three pressure taps 
are located chordwise across the upper surface and two on the lower 
surface at midspan. Figure 6 is a schematic drawing of the model showing 
the location of the pressure taps. 
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FIGURE 6. PRESSURE TAPS ON WIND TUNNEL MODEL 













The model was mounted spanwise across the test section on a mechanism 
which allowed measuring the angle of attack about the midchord line. The 
pressure tap leads were brought through the tunnel walls with stainless 
steel tubing. Figure 7 shows the model mounted in the tunnel test section. 

C. PRESSURE MEASURING SYSTEM 
1. Pressure Transducers 

Limitations on size and expense lead to the selection of a remote 
pressure transducer designed by Professor L. V. Schmidt of the Naval 
Postgraduate School. The design is shown in Figure 8. Two transducers 
were built, one to measure pressures on the model and the other to be 
used as a reference. 

The transducer system consisted basically of a Bently Detector 
System, a 0.003 inch thick annealed brass diaphragm mounted in an 
aluminum housing, and plastic and steel tubing enabling the pressure on 
the airfoil to be transmitted to the transducer diaphragm cavity. The 
diaphragm thickness was selected so as to provide a useable instrument 
in the pressure range of + 0.75 psia. The disassembled transducer is 
shown in Figure 9. 

The Bently Detector System provides a voltage signal which is 
linear with respect to the distance change detected by its probe. For 
small pressure differentials across the brass diaphragm, the deflection 
is proportional to pressure. By putting a known pressure differential 
across the diaphragm, the transducer can be calibrated yielding a 
linear calibration curve for small differentials. The details of this 
procedure are given in Appendix A. 

The stainless steel and plastic tubing act as transmission lines 
from the pressure taps to the transducer. The unsteady pressure components 
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FIGURE 8.SECTIONAL DRAWING OF PRESSURE TRANSDUCER 
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FIGURE 9. PHOTOGRAPH OF TRANSDXER COMPONENTS 
















undergo frequency-dependent phase shift and attenuation through the 
tubing. A dynamic calibration (Appendix A) was performed at 5 hertz 
increments and intermediate values at 1 hertz increments were obtained 
by linear interpolation and extrapolation. The results may be seen in 
Figures 10 and 11. 

2. Power Supplies and Amplifier 

A Power Design Transistorized Power Supply, Model 205IR was used 
to power the Bently Proximeter. The Proximeter output at zero pressure 
differential was in the neighborhood of + 5.73 Volts as installed in the 
pressure measuring system. A Kepco Model OPS 21-1 Operational Power 
Supply/Amplifier was used to amplify the Proximeter signal. The output 
of the operational amplifier was adjusted to zero volts at zero pressure 
differential by means of a suitable bias voltage. The gain of the 
operational amplifier at a given feedback resistance was verified to be 
constant with input frequency over the frequency range of interest. The 
phase shift due to the operational amplifier was a constant -180 degrees 
over the same frequency range. Figure 12 is a schematic diagram of the 
pressure measuring system and Figure 13 is a typical static calibration. 
Figure 14 shows the pressure measuring system installed. 

D. HOT WIRE ANAMOMETER 

A constant temperature, transistorized, hot wire anemometer was used 
to measure the free stream velocity. The electrical output of the hot 
wire anemometer was proportional to the velocity. A micromanometer 
connected to a pitot-static tube was used to measure the mean free 
stream velocity and to statically calibrate the hot wire anemometer. 
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FIGURE O. TRANSDUCER GAIN AS A FUNCTION OF FREQUENCY 















Frequency, Hertz 



33 


FIGURE II. TRANSDUCER PHASE ANGLE, 
AS A FUNCTION OF FREQUENCY. 
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FIGURE 12. SCHEMATIC DIAGRAM OF PRESSURE MEASURING SYSTEM. 
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FIGURE 13. TYPICAL STATIC CALIBRATION 
CURVE OF PRESSURE MEASURING SYSTEM 
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E. TAPE RECORDER 


The outputs of the hot wire anamometer and pressure transducer were 
recorded on an Ampex CP-100 Portable Instrumentation Magnetic Tape 
Recorder/Reproducer. Typical record and reproduce calibration values 
for the recorder at a tape speed of 3-3/4 inches per second are shown 
in Table I. 

TABLE I 

TYPICAL TAPE RECORDER CALIBRATION 
Record 

D.C. INPUT (VOLTS) FREQUENCY (HERTZ) 

1.414 9448 

0 6748 

-1.414 4047 


Reproduce 


FREQUENCY INPUT (HERTZ) D.C. OUTPUT (VOLTS) 

9450 1.420 

6750 0.0 

4050 -1.409 

F. MISCELIANEOUS 

The velocity and pressure analogues were monitored on a Tektronix Dual 
Beam Oscilloscope. Pictures of representative velocity and pressure 
outputs were recorded photographically from oscilliscope traces. 

The shutter valve frequency was measured with a magnetic pickup 
located outboard of the upper shutter valve blade shaft. The output of 
the pickup was read on a decade counter. A schematic diagram of the 
complete data aquisition system is shown in Figure 15. 
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FIGURE 15. SCHEMATIC DIAGRAM OF DATA 
ACQUISITION SYSTEM ‘ 
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IV. PROCEDURES 


Prior to each run, a static calibration of the transducer was made in 
order to ensure that no drift had occurred since the previous calibration. 

The desired free stream mean velocity was established in steady flow 
and the hot wire anemometer set to give the desired output at the free 
stream mean velocity. The oscillations were then adjusted to the required 
frequency. 

The angle of attack vernier was calibrated by comparing the trans¬ 
ducer output of the taps on the lower surface with those at corresponding 
locations on the upper surface. The angle of attack vernier was accurate 
to + .05 degrees when calibrated on this basis at zero degrees angle of 
attack. 

Apparent regularity of the free stream wave-form on the oscilloscope 
was used to establish the desired oscillating frequencies. Subsequent 
Fourier analysis of the free stream velocity wave-form showed this pro¬ 
cedure to be adequate in establishing clean wave-forms. 

The airfoil was placed at the desired positive angle of attack and 
recording began. The output of the transducer connected to each pressure 
tap in turn was recorded on tape simultaneously with the output of the 
hot wire anamometer. Since the airfoil was symmetrical, the pressure on 
the lower surface was assumed to correspond to the pressure on the upper 
tapped surface at negative angle of attack. The airfoil was rotated to 
a negative angle of attack in order to obtain the pressure distribution 
on the lower surface. 
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V. DATA ANALYSIS 


A. ANALOGUE TO DIGITAL CONVERSION 

The recorded data were converted from analogue to digital repre¬ 
sentation on the Naval Postgraduate School's Hybrid Computer. Details 
of the data conversion process are given in Appendix B. 

The analogue signals were filtered through two low-pass R. C. filters 
in series prior to digitization. The filter was used to improve signal 
to noise ratios of the analogue data. Preliminary investigations had 
shown that some high frequency noise was picked up in the data recording 
process from both the transducer and the tape recorder. ("High frequency", 
for the purposes of this investigation, means frequencies greater than 500 
hertz.) The filters were built around the existing high gain amplifiers 
on the analogue side of the hybrid. The filter transfer function was 
measured at 10 hertz increments and intermediate values at 1 hertz 
increments obtained by linear interpolation. Figures 16 and 17 give 
the transfer function of the filtering system. 

All sampling was done at the rate of 1024 samples per second. The 
analogous pressure and velocity signals associated with each pressure tap 
were sampled for 4 seconds for a total sample size of 4096 samples per 
digital record. 

Prior to each run a D. C. reference voltage measured with a digital 
voltmeter was recorded. These reference voltages were then digitized as 
a check on the tape recorder and digitizing process. Typical values of 
voltmeter measured D. C. voltages and the corresponding digital mean 
value are shown in Table II. 
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RGURE 17. FILTER PHASE ANGLE AS A FUNCTION OF FREQUENCY 








TABLE II 


RECORDED D.C. VOLTAGE AND CORRESPONDING DIGITAL MEAN VALUE 


RECORDED (VOLTS) 


DIGITAL MEAN VALUE (VOLTS) 


+ 1.414 


+1.409 

0.995 

0.496 

0.104 

0.002 

- 0.100 

-0.496 

-0.994 

-1.415 


1.0 

0.5 

0.1 

0.0 

- 0.1 

-0.5 

- 1.0 


-1.414 


B. TAPE CONVERSION 

The arrays of digitized data obtained from the analogue to digital 
conversion process were stored on 7 track, 200 bit per inch magnetic 
tape in 24 bit words. A conversion to 9 track, 800 bit per inch tape 
and 32 bit words was necessary in order to analyze the data on the IBM 
360 Series Computer. The computer program used to make the conversion 
is included as an appendage to this thesis under the title "TAPE CON¬ 
VERSION PROGRAM". 

C. ADDITIONAL ASSUMPTIONS 

Ideally, phase information between the free stream velocity and 
pressure analogues at a given tap would have been preserved during the 
analogue to digital conversion process. Because of equipment limitations 
discussed in Appendix B, this was not possible. However, measurements 
from photographs of oscilloscope traces of the velocity and pressure 
analogues corresponding to given pressure taps indicated that the 
phasing between the free stream velocity and the pressure on the surface 
of the airfoil was essentially independent of the surface coordinate. 

This measured time difference between a negatively-sloped zero crossing 
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of the corresponding pressure analogue was therefore taken as constant 
for a given run. 


D. NON-DIMENS lONALIZATION 

1. Pressure Coefficients 

The pressure data corresponded to the difference between atmos¬ 
pheric pressure and the surface pressure. It was therefore convenient 
to define a non-dimensional pressure coefficient referenced to atmos¬ 
pheric pressure as: 

P(x,y,t)-P 

Cp(x.y,t)-^ 

q 

where : 


P(x,y,t) = pressure at time t at a point (x,y) on the surface of 
the airfoil, 

P = atmospheric pressure, and 

A 

q = mean free stream dynamic pressure. 

The mean pressure coefficient, Cp(x,y), is then defined as: 

P(x,y)-P 

Cp(x,y) = -- and 

q 

the ’’unsteady pressure coefficient*’ is defined as: 


Cp(x,y,t) 


_ y 

q 


where : 


p(x,y,t) = "unsteady pressure’’ at time t at a point (x,y) on the 
surface of the airfoil. 

2. Normal Force Coefficients 


The time dependent force per unit span acting normal to the chord 
of a two-dimensional airfoil is approximated by: 
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c 


N(t) = J [P(x,-y,t)-P(x,y,t)] dx 
0 

if the airfoil is symmetrical. The normal force coefficient can be 

defined in terms of the non-dimensional pressure coefficient as: 

1 

Cj^(t) = J [Cp(x,-y,t)-Cp(x,y,t)] d§ 

0 


where ; 

§ = x/c. 

In terms of the mean and unsteady pressure coefficients this becomes: 

1 1 

C^(t) = /[Cp(§,-Tl) - Cp(§,in)] d^+ J [Cp(§,-Ti,t)-Cp(§,Ti,t)] d§ 

0 0 


which leads to: 




I- + 

qc qc 



-f c^(t) 


where : 


C-, = mean normal force coefficient 
N 

c^(t) = "unsteady normal force coefficient". 

3. Center of Pressure 

The center of pressure may be defined in terms of non-dimensional 
parameters in a manner similar to that used above. The result is: 


1 


I [Cp(§,-T1,t) - Cp(§,r,t)] § d§ 

F = 9 . - 

^c.p. C^(t) 

E. COMPUTATIONAL PROCEDURES 

The numerical techniques used in the calculations are discussed in 
Appendix C. A flow chart of the computer program used to analyze the 
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digitized data is shown in Fig. 18. A listing of the program is included 


as an appendage to this thesis. 
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FIGURE 18. ABBREVIATED COMPUTER PROGRAM FLOW CHART 
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VI. RESULTS AND DISCUSSION 


A. CONDITIONS INVESTIGATED 

Table III lists the freestream conditions and angles of attack at 
which data was taken. A *'Run” consisted of recording the pressure at 
each tap for 30 seconds with the airfoil at a given angle of attack in 
a free stream velocity with a given mean velocity, amplitude of oscil¬ 
lation, and frequency of oscillation. 

The steady flow Runs were made and the pressure distribution and 
normal force coefficients presented in order to provide a basis for 
comparison with the non-steady results. Power spectral densities of 
the turbulent pressure fluctuations in steady flow were calculated to 
investigate the possibility of any harmonic coupling that might occur 
between the forced oscillations of the free stream and the vortex 
shedding phenomena observed in steady flow. 

B. STEADY FLOW RESULTS 

The measured pressure distributions on the airfoil at angles of 
attack of 5, 10, and 15 degrees in a steady flow of approximately 100 
feet per second are shown in Fig. 19 through 21. Figure 22 shows the 
resulting normal force coefficients. The lift curve of a NACA 63-009 
airfoil section is included for comparison [Ref. 12]. 

Power spectral densities of the turbulent pressure fluctuations at 
representative taps are shown in Fig. 23 through 32. These power 
spectrums and subsequent power spectrums shown have been normalized with 
respect to their mean squares. 

At angles of attack of 0 and 5 degrees a distinct periodic component 
in the neighborhood of 21 to 22 hertz was discernable at all chordwise 
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locations investigated. The corresponding Strouhal Number taken with 
respect to chord and the frequency of the periodic component is: 

Nst = 0A05. 

At 10 degrees angle of attack the spectrum of the fluctuating 
pressure on the upper surface appears to be broadened over the lower 
frequencies. Periodic components are discernable along the lower sur¬ 
face at about the same frequency as was noted at lower angles of attack. 

At 15 degrees angle of attack the airfoil was fully stalled as 
evidenced by the pressure distribution. Figures 30 through 32 show a 
distinct periodic component between 21 and 22 hertz just as the runs at 
lower angle of attack. This periodicity was visible at all pressure 
taps on the surface. The tendency of the spectrum to spread more 
towards lower frequencies at downstream stations on the upper surface 
is apparent when Fig. 30 and 31 are compared. 

Figures 33 and 34 show the mean square of the fluctuating pressure 
coefficients for Runs 3 and 4 plotted as a function of chord. Discussion 
of the possible significance of the resulting distribution will be 
deferred to the unsteady flow section. 

Noise at 60 hertz was filtered from the spectral and mean square 
distributions shown in Fig. 23 through 34 by setting the computed value 
of the spectrum at 60 hertz equal to the average of the values at 59 and 
61 hertz. A typical unfiltered spectrum is shown in Fig. 35. 

C. UNSTEADY FLOW RESULTS 

1. Run 5 

Run 5 was made at a fundamental frequency of 94 hertz, an angle 
of attack of 15 degrees, and a free stream velocity with mean of approxi¬ 
mately 100 feet per second and amplitude ratio of 0.085. Figures 36 and 
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37 show typical oscilloscope traces of the velocity and pressure analogue 
wave-forms at selected pressure taps for this Run. 

It is apparent from these photographs that the time between a 
negatively sloped zero crossing on the velocity analogue wave and a 
positive sloped zero crossing on the pressure analogue wave is essentially 
constant at a given tap even though the wave-form is obviously not the 
same for every period. These same pictures show that the change in the 
above time-difference with chordwise location is small. 

The time-difference mentioned above was measured at two sample 
time-increments from photographs similar to those shown in Fig. 36 and 
37, but on an expanded time scale. Expressed in terms of a phase angle 
at the fundamental frequency of 94 hertz, the two sample time-increments 
are equivalent to 66 degrees. This measurement served to establish the 
phasing that was used in reducing the data. 

Figure 38 shows the power spectrum of the oscillating free 
stream velocity for Run 5. The spectral quality of the oscillations 
produced by the tunnel is immediately apparent from the manner in which 
the fundamental frequency dominates the spectrum. 

Figures 39 through 64 are power spectral densities of the unsteady 
pressures associated with the indicated chordwise locations. They were 
selected to indicate the complete trend of the pressure spectrum on the 
surface of the airfoil. The spectrums appear to be essentially the same 
along the upper surface to about 30 per cent chord at which point the 
spectral component at the fundamental harmonic begins to decrease and the 
second and third order harmonic components increase. These effects appear 
to be maximized at or near the mid-chord point where the second order 
harmonic dominates the unsteady pressure. (It is unfortunate that the 
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spacing between pressure taps near the mid-chord point increases but the 
observed pressure spectrums were not anticipated before the model was 
constructed.) Less than 20 per cent chord downstream from the mid-chord 
point the spectrum again takes on the appearance noted near the leading 
edge . 

Figure 65 shows the measured mean square of the unsteady pressure 
coefficient at each tap. There is a decrease in the mean square as mid¬ 
chord is approached and a definite increase as the downstream edge is 
approached, implying some sort of spatial correlation between the mean 
square of the pressure fluctuation and chordwise distance downstream. 

Figure 66 shows the calculated time dependent normal force 
coefficient and the free stream velocity plotted with common time origins 
as determined by the previously described photographic measurement. 

Figure 67 is the power spectrum of the calculated unsteady normal 

force. 

Figure 68 shows the mean pressure distribution as calculated for 
Run 5. Figures 69 and 70 are photographs of the time dependent pressure 
distribution displayed on an IBM 2250 Visual Display Unit. The computer 
program used to create the display is attached as an appendage under the 
title ’’DISPIAY PROGRAM”. 

Figure 71 is the power spectruia of the calculated unsteady 
center of pressure for Run 5. 

2. Run 7 

Run 7 was made at a fundamental frequency of 11 hertz, an angle 
of attack of 15 degrees, and a free stream velocity with a mean of approxi¬ 
mately 100 feet per second and amplitude ratio of 0.072. Figures 72 and 
73 show typical oscilloscope traces of the velocity and pressure analogue 
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wave-forms at selected pressure taps for Run 7. Note the regularity of 
the pressure analogue wave-forms with time for a given tap and the 
similarities of the wave-forms for the different taps. The time dif¬ 
ference between a negatively sloped zero crossing on the velocity 
analogue wave-form and a positively sloped zero crossing on the pressure 
analogue wave-form was measured at approximately 24 sample time-increments. 
At the fundamental frequency of 11 hertz, 24 sample time-increments is 
approximately given in terms of a phase angle as 93 degrees. This measure¬ 
ment was used as described previously to establish the phasing between the 
free stream velocity and pressure in data reduction. 

The power spectral density of the unsteady free stream velocity 
for Run 7 is shown in Fig. 74. The quality of the oscillation should 
again be noted. 

Power spectrums of the unsteady pressure associated with repre¬ 
sentative taps are shown in Fig. 75 through 80. The wave-form similarity 
noted in the oscilloscope pictures is also apparent in the power spectrums. 
The similarity of the spectrums to each other should also be noted. 

Figure 81 gives the mean square of the unsteady pressure coefficient 
associated with each tap as a function of the non-dimensional chord. 

The unsteady velocity and normal force coefficient are given as 
functions of time from a common time origin in Fig. 82 and 83. The power 
spectrum of the unsteady normal force is given in Fig. 84 and the power 
spectrum of the center of pressure in Fig. 85. Figure 86 shows the mean 
pressure distribution on the airfoil for Run 7. 

D. COMPARISON OF RESULTS AT DIFFERENT FREQUENCIES 
1. Normal Force and Center of Pressure 

The most obvious difference between the results presented is in 
the mean normal force coefficients. This is illustrated in Fig. 87 which 
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shows the mean normal force coefficient at 15 degrees angle of attack as 
a function of frequency. Recall that the mean velocities for the steady 
flow Run and the oscillating Runs 5 and 7 were all about 100 feet per 
second. Run 5 had a fundamental frequency of 94 hertz while Run 7 had 
a fundamental frequency of 11 hertz. The amplitude ratio for the two 
oscillatory Runs was about the same, that for Run 5 being 0.085 and that 
for Run 7 being 0.072. 

The harmonic composition of the unsteady normal force coefficient 
for the two oscillatory Runs under consideration was also much different. 

A comparison of Fig. 67 and Fig. 84 shows that the second harmonic of 94 
hertz is much more predominant than the second harmonic of 11 hertz. 

Center of pressure calculations also reflect the difference in 
harmonic composition that was found in the normal force results. The 
unsteady center of pressure at 94 hertz is dominated by the second harmonic 
and the first and third harmonics are about the same, while at 11 hertz 
little difference between the spectrum of the unsteady normal force and 
unsteady center of pressure was noted. 

2. Pressure Distributions 

The spectral composition of the unsteady pressure reflects the 
difference noted in the normal force and center of pressure results. The 
presence of higher order harmonics is apparently not just related to 
increasing frequency. Allen (Ref. 11] showed oscilloscope traces of the 
pressure analogues associated with the same airfoil used in this investi¬ 
gation at a fundamental frequency of 128 hertz and about the same mean 
and mean square velocity. His traces did not give an indication of any 
significant higher order harmonics. 

There are apparent spatial periodicities evident in the instan¬ 
taneous pressure distributions shown in Fig. 69 and 70 for Run 5. The 
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mean pressure distribution for Run 7 shown in Fig. 86 also clearly 
indicates a spatial periodicity. The mean square distributions of the 
unsteady pressure plotted in Fig. 65 and 81 show a similar trend. 

A hint as to the cause of these periodic trends may be obtained 
by considering the chordwise distribution of the mean square of the 
pressure fluctuation in steady flow. One would normally expect vortices 
shed at or near the separation point on the airfoil to immediately begin 
to dissipate as they travel downstream from their point of origin. A 
resulting decrease in the mean square of the fluctuations in velocity 
and pressure would be expected with chordwise distance downstream as is 
shown in Fig. 33 for the airfoil at 10 degrees angle of attack in steady 
flow. The mean square of the pressure fluctuation on the airfoil at 15 
degrees angle of attack in steady flow shows a definite increase toward 
the trailing edge as shown in Fig. 34. This suggests that the spatial 
periodicity observed in the unsteady runs was a result of the combination 
of the oscillatory effects and the effects of the stalled airfoil. 

The instantaneous pressure distributions shown for Run 5 in Fig. 

69 and 70 also suggest that the second harmonic component of the normal 
force could be due to a delay of separation near the leading edge associ- 
ated with the minimum of the velocity or change in sign of the acceleration 
of the free stream. Reference to the instantaneous pressure distribution 
associated with the minimum of the velocity demonstrates the more favorable 
spatial pressure distribution that occurs at that point on the velocity 
wave form. 

E. COMPARISON OF UNSTEADY FLOW RESULTS WITH STEADY FLOW ANALYSIS 

Steady flow results have traditionally been used to estimate aero¬ 
dynamic forces acting on airfoils in unsteady environments. Such analysis 
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usually involves assuming that the lift coefficient of an airfoil section 
is constant at a given angle of attack within a given range of Reynolds 
Number. 

The normal force coefficient found for the airfoil used in this 
investigation at 15 degrees angle of attack and a Reynolds Number based 
on chord of approximately 3 x 10^ was 0.681. Assuming the normal force 
coefficient to be constant for a given angle of attack leads to the fol¬ 
lowing estimates for the mean normal force coefficient and mean square of 
the unsteady normal force coefficient for the same airfoil at 15 degrees 
angle of attack in an oscillating flow with amplitude ratio e = 0.08; 

- 0.683 
s 

s 

where is the normal force coefficient in steady flow, 

s 

Note that assuming a constant normal force coefficient leads to 
results that depend only on the amplitude ratio of the oscillating 
free stream and are independent of frequency. 

The experimentally measured mean normal force coefficients for the 
airfoil at 15 degrees angle of attack in a free stream oscillating at 11 
and 94 hertz were 0.851 and 1.401 respectively. The mean squares of the 
unsteady normal force coefficients for the two cases were 0.075 and 0.100 
respectively. The difference between the experimentally derived values 
and the values predicted by assuming a constant normal force coefficient 
at a given angle of attack are considered by the author to be significant. 

F. SUMMARY OF RESULTS 

Figure 87 gives the measured mean normal force coefficient as a 

function of frequency for the airfoil at an angle of attack of 15 degrees 
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in a free stream with a mean velocity of approximately 100 feet per 
second and amplitude ratio of approximately 0.08. Figure 88 gives the 
measured mean center of pressure under the same conditions. 

Due to lack of computing time and a machine malfunction in the 
digitization process, the data taken for Runs 6 and 8 were not reduced. 
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VII. CONCLUSIONS 


From the results presented, the following conclusions may be drawn 
for the conditions investigated: 

1. The magnitude of the pressure force acting normal to the 
chord of an airfoil at high angle of attack in oscillating flow is 
frequency dependent and is significantly greater than that observed in 
steady flow and is not adequately predicted by quasi-steady analysis. 

2. Higher order harmonic components of the fundamental free 
stream frequency of oscillation are of significant order of magnitude in 
the normal force and their relative magnitude is frequency dependent. 

3. Distinct periodic components are observable in the pressure 
acting on an airfoil in steady flow. 
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TABLE III 


EXPERIMENTAL CONDITIONS FOR WHICH DATA WERE TAKEN 


Run No. 

Angle 

of 

Attack 

(degrees) 

Mean 

Ve locity 
(ft/sec) 

Frequency 

of 

Oscillation 
(hertz) 

Amplitude 

Ratio 

of 

Oscillation 

1 

0 

100 

steady 

-- 

2 

5 

100 

steady 

-- 

3 

10 

100 

steady 

-- 

A 

15 

100 

steady 

-- 

5 

15 

100 

94 

0.08 

6 

15 

150 

94 

0.08 

7 

15 

100 

11 

0.08 

8 

15 

70 

9 

0.25 
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FIGURE 25. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 2 
^'00 MEAN SQUARE =0.00334 
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Frequency, Hertz 

FIGURE 26. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 2 
(.-OA5 UPPER SURFACE MEAN SQUARE=0.00993 
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FIGURE 27. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 3 
^=0.0 MEAN SQUARE -0.0727 
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FIGURE 29. PSD OF UNSTEADY PRESSURE COEFRCIENT, RUN 3 
<^ = 0.5 LOWER SURFACE MEAN SQUARE =00II6 
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^=0.125 LOWER SURFACE MEAN SQUARE = 0.0104 
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FIGURE 33. MEAN SQUARE OF UNSTEADY PRESSURE 

COEFFICIENT ON AIRFOIL 

RUN 3 (X = 10° V«lOOft / sec STEADY FLOW 
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FIGURE 34. MEAN SQUARE OF UNSTEADY PRESSURE 
COEFFICIENTS ON AIRFOIL 

RUN4 0 = 15® U«IOOft/sec STEADY FlOW 
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Figure 36. Selected Oscilloscope Traces, Run 5 
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Figure 37. Selected Oscilloscope Traces, Run 5 


77 


























ii 

0) 

X 


>\ 



00 

fO 

U- 


78 


PSD OF UNSTEADir VELOCTTY, RUN 5 
U = I02.66 V*38.43 
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FIGURE 39. PSD OF UNSTEADY PRESSURE CXDEFFICIENT, RUN 5 
6 =0.0 UPPER SURFACE 
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Frequency, Hertz 

FIGURE40. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
6^=0025 UPPER SURFACE 
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FIGURE 41. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
^0.075 UPPER SURFACE 
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FIGURE 43. PSD OF UNSTEADY PRESSURE COEFFICIENT ,RUN 5 
^=0.175 upper surface 
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Frequency, Hertz 

FIGURE 44. PSD OF UNSTEADY PRESSURE COEFFICIENT RUN 5 
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FIGURE 46 PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
6=0.350 UPPER SURFACE 










40 ^ 



-8 



4 



ajonbs uoai^ luaajad 


87 


FIGURE 47. PSD OF UNSTEADY PRESSURE GOEFFICIENT, RUN 5 
^=0-400 UPPER SURFACE 
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FIGURE 48, PSD OF UNSTEADY PRESSURE COEFRCIENT.RUN 5 
^=0.450 UPPER SURFACE 
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6^0.500 UPPER SURFACE 
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FIGURE 50. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
6«a600 UPPER SURFACE 
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Frequency .Hertz 

FIGURE 51. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
^ = 0.700 UPPER SURFACE 
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RGURE 52. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
t=0.750 UPPER SURFACE 
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Frequency, Hertz 

FIGURE 53. PSD CF UNSTEADY PRESSURE COEFFCIENT, RUN 5 
t.=0.0 LOWER SURFACE 
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FIGURE 54. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
^=0,050 LOWER SURFACE 
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FIGURE 55. PSD OF UNSTEADY PRESSURE COEFFICENT, RUN 5 
^ = 0.100 LOWER SURFACE 
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FIGURE 56. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
6=0.150 LOWER SURFACE 
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FIGURE 59. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 

(S = 0.350 LOWER SURFACE 
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FIGURE 58. PSD OF UNSTEADY PRESSURE COEFFICIENT . RUN 
£=0.250 LOWER SURFACE 
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FIGURE 59. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 

<S= 0.350 LOWER SURFACE 
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FIGURE 60. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN ! 
^=0.450 lower SURF/CE 
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FIGURE 61. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
i=0.600 LOWER SURFACE 
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FIGURE 62. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
^=0750 LOWE R SURFACE 













a 




8 

CVJ 




§ 


o 

iO 


ii 


Q> 

X 

& 


c 

0 ) 

3 



103 


FIGURE 63. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
^=0.850 LOWER SURFACE 
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FIGURE 64. PSD OF UNSTEADY PRESSURE COEFFICIENT. RUN 5 
^^0.950 LOWER SURFACE 
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FIGURE 67. PSD OF UNSTEADY NORMAL FORCE COEFFICIENT, RUN 5 

MEAN SQUARE =0.100 
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FIGURE 68. MEAN PRESSURE aSTRIBUnON ON AIRFaL 
<x=l5® fo = 94Hz S=a085 
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Figure 69. Time Dependent Pres.sure, Run 5 
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Figure 70 













4CH 




N 



O 

-22 



j 




cc 

i_ ro 

S8 


O 

LJ 


LlJ 

a: 

3 

CO 

CO 

LJ 

?E 


o 

II 

LJ 

q: 

< 


CO 

z 

< 

LJ 


ID 



J— 

CO 

z 

3 

u. 

o 

o 

CO 

CL 





LJ 

cr 

3 

Li. 


9Jonbs uo9^ luaojacd 


111 










UPPER SURFACE 



TAP 7 




TAP 5 


TAP 9 


TAP I 1 



TAP 13 TAP 14 
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Figure 72. Pictures, Run 7, Pressure, Velocity 
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Figure 73. Pictures, Run 7, Pressure, Velocity 
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FIGURE 74 . OF UNSTEADY VELOCITY, RUN 7 
U = 100.46 U2=25.99 
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FIGURE 75 PSD OF UNSTEADY PRESSURE COEFFICIENT, RL 
t = 00 UPPER SURFACE 
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FIGURE 76. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 7 
t = OIOO UPPER SURFACE 
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FIGURE 77. PSD OF UNSTEADY PRESSURE COEFFICIENT . RUN 7 
(f =0.500 UPPER SURFACE 
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FIGURE 79 PSD OF UNSTEADY PRESSURE COEFROENT^ RUN? 
i=0.100 LOWER SURFACE 
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Frequency .Hertz 

FIGURE 80. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 

LOWER SURFACE 
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FIGURE 82. TIME DEPENDENT VELOCrTT, RUN 7 
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NORMAL FORCE , RUN 7 
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MEAN SC3UARE =0.075 
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MEAN SQUARE =0.010 
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FIGURE-86. MEAN PRESSURE DISTRBLTTDN ON AIRFOL 
0^=15® fo=ll Hz ^=0.072 









FIGURES?. MEAN NORMAL FORCE COEFFICIENT VS. FREQUENCY 
0^=15° U=IOOft/tec £,^0.08 



FIGURE 88. MEAN CENTER OF PRESSURE VS-FREQUENCY 
0^=15° U =100 ft/sec ^5^0.08 
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APPENDIX A 


TRANSDUCER CALIBRATION 

Dyanmic calibration was carried out using the procedure developed 
by Johnson in Reference 14. The results of his work were used to select 
the inside diameter of the tubing and tube length for this investigation. 
The stainless steel tubing used in the model was 0.0625 inches in outside 
diameter, 0.047 inches in inside diameter and 24 inches long. The plastic 
tubing attached to the transducer housing was 0.49 inches in inside dia¬ 
meter and 11.5 inches long. 

Both transducers were first mounted in a calibration chamber and 
statically calibrated. Figure 89 shows a schematic diagram of the static 
calibration instrumentation. 

A sound driver was then mounted on the calibration chamber and driven 
at known frequencies by the amplified output of an audio oscillator. The 
spectrum from 0 to 500 hertz was scanned to see if any difference in phase 
or gain could be detected between the two transducers. No difference in 
phasing was detected and the ratio of the gains was constant for the 
observed spectrum. 

The next step was to remove the transducer to be used for the investi¬ 
gation from the chamber, install it in its housing and reconnect it to the 
chamber through steel and plastic tubing identical in length and diameter 
to that to be used in the experiment. The values of phase difference 
and wave analyzer readings for each transducer were then recorded at 5 
hertz increments from 0 to 400 hertz. The gain of the remote transducer 
was then taken as the product of the ratio of static gains and the ratio 
of the wave analyzer readings. Hie amplifiers used to amplify the two 


128 




transducer outputs were switched and the procedure was repeated. No 
differences in results were noted. Figure 90 shows a schematic of the 
dynamic calibration instrumentation. 

The spectral interval subsequently used in the analysis of results 
was 1 hertz. A linear interpolation was used to obtain intermediate 
points in the frequency range 0 to 400 hertz and a simple linear extrapo¬ 
lation made from 400 to 512 hertz. The results are shown in Fig. 10 and 
11 . 
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D.C. VOLTMETER 


R6URE 89. STATIC CALIBRATION INSTRUMENTATION 
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R6URE 90. DYNAMIC CALIBRATION INSTRJMENWION 















































APPENDIX B 


ANALOGUE TO DIGITAL CONVERSION 

The available digitizing equipment was not capable of digitizing 
more than one analog signal at a time. This precluded simultaneous 
digitization of the velocity and pressure analogues associated with a 
given tap and ultimately required photographic measurement of phasing. 

The procedure used was based on an available digital computer pro¬ 
gram. The program takes binary words from one channel of the analogue to 
digital interface of the hybrid and stores them in a given buffer until 
the buffer is full. The contents of that buffer are then written on 
magnetic tape while a second buffer is being filled. When the second 
buffer is full it is written on magnetic tape and the first buffer takes 
on the storage function. The cycle repeats itself for the sample length 
desired. Sample length is expressed in terms of the number of blocks of 
1024 samples per record. The 1024 samples correspond to the storage 
capacity of the buffers. 

The analogue signal to be digitized is taken by a designated trunk 
line fromthe analogue side of the hybrid to the analog to digital con¬ 
verter. The sampling interval is controlled by interupts into a trunk 
line on the logic board of the analogue. 

For this investigation a sampling frequency of 1024 samples per 
second was chosen in order to simplify the computer analysis of the 
data. For minimum scale error in the conversion process it was desirable 
that the digitized signal be as large as possible. A compromise between 
the desired amplitude of the digitized signal and numerical convenience 
was worked out by making the filter gain be of the order of 100. 
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The filter was used to improve signal to noise-plus-signa1 ratios of 
the analogue data. Preliminary investigations had shown that some high 
frequency noise was picked up in the data recording process from both 
the transducer and the tape recorder. ('*High frequency** for the purpose 
of this investigation meaning frequencies greater than 500 hertz.) 
Filtering prevented this high frequency noise from distorting the 
resulting spectral estimates. 

For the steady flow runs it was not necessary to preserve any phase 
informat ionand the analogue circuitry used is shown schematically in Fig. 

91. The conversion process in this case was commenced when the **Master 
Clock**on the analogue computer was switched from **Reset** to **Run*' with the 
analogue computer in the ^Compute** mode. In more specific terms, switching 
the master clock to **Run** enabled the interrupts to the digital computer. 

The unsteady flow runs were digitized using the circuitry shown in 
Fig. 92. A 1000 hertz sine wave was recorded on a third track of the 
data tape. This signal was used to start the interrupts to the digital 
computer at the point on the tape corresponding to the start of the 1000 
hertz signal. The idea was that the velocity and pressure analogues corres¬ 
ponding to a given pressure tap would be digitized starting at the same 
point on the tape. The large gain on the 1000 hertz signal was used to 
get maximum rise time to the comparator. Utilizing the circuitry shown 
in Fig. 92 also allowed control of the digitizing process from the 
digital side of the hybrid. 

The procedure described above was not adequate to preserve the 
phase information between the velocity and pressure analogue signals, 
Possible reasons include variabilities in the response times of the 
comparator and/or the delay flop in the logic circuitry and/or recording 
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data at too slow a tape speed to be able to get the rise time necessary 
to insure adequate responses of the above analogue devices. 

A method that was not tried would be to use the 1000 hertz signal 
to put the analogue into **Compute’* rather than utilizing the delay flop 
and comparator circuitry shown in Fig. 92. This would eliminate at 
least one of the elements in the circuit and possibly both. 

The hybrid has the capability of converting32 analogue channels with 
only a slight time skew. The conversion process is limited by buffer 
size and the speed at which the filled buffers can be written on magnetic 
tape. The time it takes to fill a buffer cannot be less than the time it 
takes to write the contents of a buffer on magnetic tape. This limits the 
sampling frequency. 

In order to sample more than one track at a time the digital computer 
program would have had to be written to store the sampled analogue signals 
alternately in the buffers. The author did not have the time available 
to rewrite the program and the Electrical Engineering Department Computer 
Laboratory does not have the programming support available to rewrite it. 
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CIRCUITRY C ANALOGUE IN COMPUTE] 
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FIGURE 92. IDGIC AND ANALOGUE CIRCUITRY FOR UNSTEADY 

FLOW CONVERSION 






























APPENDIX C 


NUMERICAL ANALYSIS OF DATA 


A. ASSUMPTIONS 

The digital analysis was based on the assumption that the pressure 
at any point on the surface of the airfoil was a parametric function of 
the free stream velocity, i.e., 

P = P(x,y,U), U = U(t) 

The oscilloscope traces shown in Fig. 36 and Fig. 72 indicate the validity 
of this assumption which allowed the establishment of a common time base 
for each pressure analogue sample. 


B. MEAN AND MEAN SQUARE CALCULATIONS 
The mean of a digital array 
relation: 




N was calculated from the 


V = 


N 

E V 
n=l ^ 


where corresponds physically to a quantity measured at time t^ and N 
is the number of samples in the array. For spectral calculations, it was 
desirable that the arrays subjected to spectral analysis have a zero mean. 
This prevents distortion at the low frequency end of the spectrum [Ref. 13]. 
Consequently, the mean was subtracted from the initial array <V > to form 
a new digital array This new array corresponded to the unsteady 

part of the analogue signal. 

The mean square of calculated from the relation: 

N 


?. i i; u ^ 

N T n 
n=l 
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C. SPECTRAL CALCULATIONS 


1. Fourier Transforms 

All sampling for this investigation was done at the rate of 1024 
samples per second. Spectral analysis was performed on arrays of 1024 
samples each. This procedure yields a Nyquist Cutoff Frequency of 512 
hertz and an equivalent spectral bandwidth of 1 hertz [Ref. 13]. 

The complex Fourier Transform of the digital array under con¬ 
sideration was calculated using the relation: 

N-1 9 . nK 

2 U . . N 

K N ^ nH-1 
n=0 


where : 

K-Fq 

= spectral component at frequency , K = 1, ... N/2 

i = V-T 
e = 2.7183... 

The inverse transfer functions of the tubing and filter were 
applied to the pressure analogues in complex form in order to get the 
Fourier Transform of the pressure at the surface of the airfoil, i.e., 

p 

«k“k 

The Fourier Transform of the unsteady velocity was found by 
simply dividing the Fourier Transform of the velocity analogue by the 
transfer function of the filter. 

The inverses of the pressure and velocity transforms were taken 
using the same subroutine as was used to take the transform of the 
original analogue array. Details of the capabilities of this sub¬ 
routine, called "HARM” in the data analysis program, are given by 
Cooley and Tukey in Reference 15. 
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2. Power Spectral Densities 



The power spectral density of 
are related by: 






and the Fourier Transform of 


The power spectral density of a function describes the general frequency 
composition of the function in terms of the spectral density of its mean 
square. Thus the mean square of was calculated from: 

00 




= Z G, 
K=0 


Af. 


K 


where the equivalent spectral bandwidth, Af , was 1 hertz for this 

D 

investigation. Reference 13 has a complete discussion of the properties 
of the power spectral density and its uses. 


D. NORMAL FORCE AND CENTER OF PRESSURE 

The integral relations given in Chapter IV, D, were approximated by 
a simple trapezoidal summation. For example, the mean normal force 
coefficient, defined as 

1 

Cn =/ [Cp (%,-r]) - Cp (§,+ 1 -)] d§ 

0 


was approximated by: 


C 


N 


23 (AC (^1) + AC,(n)l 
^ ^- 2 - 


where: 


ACp(n) = Cp (n) - Cp (n). 

u L 
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APPENDIX D 


EVALUATION OF EXPERIMENTAL TECHNIQUES 

A. PRESSURE MEASURING SYSTEM 

Dynamic calibration of the transducer and pressure tap leads was 
accomplished using a sound driver to produce reference pressure fluctu¬ 
ations, (Appendix A). The lower frequency bound on the flat portion of 
the dynamic response curve of the sound driver was 160 hertz with 
detectable response of decreasing amplitude to 25 hertz. 

Below 25 hertz, both the gain of the amplifier used to drive the 
sound driver and the gain of the sound driver were very small. The 
power required to produce any output response from the sound driver 
below 25 hertz caused saturation of the amplifier with the result that 
the output of the sound driver was a distorted signal rather than a 
simple sine wave. The combination of a complex signal, decreasing 
interval between harmonics, and low gain made the wave analyzer more 
difficult to read and the resulting gains subject to greater uncertainty 
than gains measured at higher frequencies. 

Defining error in terms of a maximum percentage of deviation from a 
true value, the estimated error in the system transfer function was 5 
per cent above 25 hertz and 10 per cent below 25 hertz. 

B. TAPE RECORDER AND ANALOGUE TO DIGITAL CONVERSION 

Table II gives examples of the differences between measured recorded 
voltages and their digital equivalents. The differences were considered 
to be insignificant. 
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C. SPECTRAL ANALYSIS 


The power spectral density estimates of the unsteady pressures and 
velocities are considered by the author to be reasonable estimates to 
their true power spectral densities, A detailed discussion of error 
criteria in spectral calculations is given by Bendat and Piersol in 
Reference 13. 

D. NORMAL FORCE AND CENTER OF PRESSURE 

The inability to preserve phase information in the analogue to 
digital conversion process and subsequent measurement of time differ¬ 
ences from oscilloscope traces made error in the normal force calculations 
difficult to evaluate. 

It was estimated that phasing between the pressure and velocity 
analogues associated with a given tap could be measured to within 1/2 
sample time increment or approximately 17 degrees at 94 hertz and 2 
degrees at 11 hertz. The quality of these measurements would establish 
minimum error bounds on the normal force and center of pressure calculations. 

E. EXPERIMENTAL ENVIRONMENT 

The variation in the ratio of the unsteady to steady components of 
the velocity associated with individual taps for Runs 5 and 7 are given 
in Tables IV and V. The pressure coefficients associated with a given 
tap were calculated using the measured mean velocity associated with that 
tap in an effort to minimize the effect of the variations on the normal 
force and center of pressure calculations. These variations in free 
stream velocity apparently were due to fluctuations in line voltage and 
atmospheric disturbances. 
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TABLE -TM. 

RATIO OF SORT(2)*SQRT(VELOCITY MEAN SQUARE) TO MEAN VELOCITY 


r< uiN 
TAP 

LOCATION 

UPPER 

LOWER 

1 

O.C 

0. ">88 

r .078 

2 

0.025 

C.091 

r.081 

3 

0.050 

0.090 

0.080 

4 

0.075 

0.089 

r.C80 

5 

0. 100 

0.089 

0.080 

6 

0.125 

0.089 

0 .0 80 

7 

0.150 

0.086 

C.083 

fi 

0.175 

0.090 

C,C83 

9 

0.200 

0.088 

0.080 

ir 

0.225 

0,090 

0,079 

11 

C.250 

0.088 

0.080 

12 

0.30C 

C.^87 

0,081 

13 

0.350 

0.091 

0 .C81 

14 

0.400 

0.088 

0.081 

15 

0.450 

0.089 

0 .082 

16 

C.500 

0.091 

0.0 82 

17 

0.600 

0.087 

0 .C81 

18 

0. 700 

O.ooi 

C.084 


0.750 

0.090 

0.086 

20 

0. R-^O 

0.087 

0.0 94 

21 

C. 850 

0.086 

0.086 

22 

0,900 

0.0R9 

0.0 86 

23 

C.O50 

0.087 

0.087 

24 

C.250 

0.088 

0.086 

25 

0.500 

C.090 

0,086 

THE 

AVERAGE MEAN 

VELOCITY FOR RUN 5 

WAS 

102,66 



THE 

AVERAGE MEAN 

SQUARE WAS 

38.43 

THE 

RATIO OF UNSTEADY TO STEADY WAS 0 

THE 

AVERAGE Q WAS 

1 2 .5 
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TABLE If 

RATIO OF SORT(2)*S0RT(VELOCITY MFAN SQUARE) TO mFAN VELOCITY 
RUN 7 


1 

0.^ 

0.069 

C .073 

2 

0.025 

0. 072 

0.071 

3 

0.050 

0.072 

0 .0 73 

4 

0.075 

0.071 

0 .072 

5 

0. IOC 

0.070 

C . 0 7 3 

6 

0.125 

0.072 

0 .0 71 

7 

0.150 

0.072 

C . C 71 

8 

0.175 

0.074 

0.0 72 

9 

0.200 

0.071 

0.073 

1C 

0.225 

0.072 

0.071 

11 

0.250 

o..:7i 

0.071 

12 

C.300 

0.070 

0 .071 

13 

0.350 

0.071 

0.074 

14 

0.400 

0.073 

^ .f 71 

15 

0.450 

0.071 

" .073 

16 

0.500 

0.071 

6.671 

17 

0.600 

0.6 74 

0 .071 

18 

0.700 

0.072 

0.071 


0.750 

0.069 

0.0 71 

2r 

0. 800 

C.070 

0 .0 7 2 

21 

0.850 

0.073 

O.C73 

22 

O.OOO 

0.071 

0.6 74 

23 

C. 950 

0.070 

0.v^7? 

24 

0.250 

0.071 

0.072 

25 

C.50C 

0.0 72 

C .074 


THE AVERAGE MEAN VELOCITY FOR RUN 
WAS 1C'!.46 

THE AVERAGE MEAN SQUARE WAS 25.99 
THF RATIO OF UNSTEADY TO STEADY WAS 
THE AVERAGE Q WAS 11.63 


0.0718 
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EXPERIVFNTAL INVESTIGATION 
INTO THE 

UNSTEADY NORMAL FORCE 
ON AN 

AIRFOIL IN AN OSCILLATING FREE STREAM 
DRCGRAMMEC BY M. R. BANNINGt CART. IJSMC. AUG-SEP. 1969 


C 


C 

C 


C 


C 


C 

C 


DI MENS ION P{4C96) ,U(6.C 96) ,CN( 10 24 ), BMP (2,2 5) , 
*BMSP(2,25) ,3MU(2,25),BMSU(2,25) ,DATA(1024) , S(256 ) , 
■!‘INV(256),M(3),Y(5J).ETA(25),GX(512),X(25) , KPHA SE ( 2 t 2 5 ) 
*FPFQ(512 I,TIMEIT(512) ,RANGE(4) 

COMPLEX’! R A( 1024, 1 , 1) ,G{ 512) ,H( 512 ) ,CON 
INTEGER SIOE,HOLE,SF,rRUN,ALPHA 


ITIL5(12)/ 


REAL’^8 ITIL1(12)/ 
^RUN t 

Rf AL-B ITIL2( 12)/ 
=^PUN 6 

REAt=^fc ITIL3(12)/ 
••!'C'UN 6 

REAL*8 ITIL4{12)/ 
’■’RUN 6 

K) AL’^e 
=«-RUN 6. 

RCAL»:=b ITIL6(1?)/ 
’I’RGN 6 

Ri:AL=:<8 1 TILTH 2)/ 
’^RUN 6 

PEAL>:=P ITIL:3(12)/ 
*RUN E 

PEAL LABEL /4H 

REAL^fl LABEL /8H 

REAL L ABC (53) 


•UNSTEADY PRESSURE COEFFICIENTS VS TI 

M.R.BANNING*/ 

•UNSTEADY VELOCITY VS TIME 

M.R.BANNING*/ 

•PSD OF UNSTEADY PRESSURE 

M.R.BANNING*/ 

•PSD OF UNSTEADY VELOCITY 

M.P.BANNING*/ 

•MEAN PRESSURE COEFFICIENTS VS ETA 

M.P.BANNING*/ 

•PSD OF UNSTEADY NORMAL FORCE 

M.R.BANNING*/ 

•UNSTEADY NORMAL FORCE VS TIME 

M.R.BANNING*/ 

•PSD OF UNSTEiDY CENTER OF PRESSURE 

M.R.BANNING*/ 

/ 

/ 


CFFINF FILE 1( 51200,2,U, IB ) 

READ(5,5(' 1 ) NRUN,NBSKIP, IF F IL , FP UN , AL PHA 
PFAC(5,503) PSLOPE,UMVQLT,QFO,TFMP,PATM 
READ(5,502){FTA(I) ,1 = 1,25) 

PEAn(5,5C4) {G( I), 1=1,512) 

Rr.AD(5,5r4) { HI I), 1 = 1,512) 

READ(5,505)(LABC( I ) , I = 1,50) 

GASCCN = 2116.0/ ( ^ 23 7 6 9-''5 1 8.6 8 8 ) 

TEMPR=TEMP+45R.688 
RH0=PATM*7C.73/(GASCnN*TEMPR) 

UMEAN=SORT(2.0*GF0*5.202/(2.54*rhu)) 


RPITF(6,7r1) NRUN,NBSKIP,PSLUPE,UMVnLT,UMEAN 
701 FCRMATIIHI, ///, T32 , • INPUT PARAMETERS FOR RUN*»I3,/, 
’^T32,*NO. OF BLOCKS SKIPPFD BEFORE READING FIRST*, 

’^* PFCCRO=* , 14,/, 

^T32,'TRANSDUCER SLCPE=•,F5.3,* VOLTS/IN. OF H20',/, 
*T32,'HnT RIRE OUTPUT OF •,F3.1,* VOLTS CORRESPONDS TO* 
•'/,T35,*A VFLCCITY OF *,F6.2,* FFET/SEC) 

KPITFI6,7C 2) FRUN,TEMP,PATM,RHC,UMEAN 
7C2 F( RMAT(T32,* FREE STREAM OSCILLAT IUN=•,14,• HERTZ*,/, 
’!'T32 ,'A TMOSPHER IC T EMPE RATURE= • , F5 . 1, • DEGREES F*,/, 
>T?2, 'ATMOSPHER IC PRESSURE=•,F6.2,' IN. OF HG*,/, 

*T32,*ATMOSPHER IC DENS ITY=•,F3.5,• SLUGS/FT',/, 
*T32,'MFAN VELOCITY=*,F7.2,• FT/SEC) 

WPITF(6,703) ALPHA 

703 FCPMAT(T32,*ANGLE CF ATTACK=*,I?,* DEGREES') 


RtWIND 4 


NCBLKS=‘^‘ 

KIN=14C 
I rSET = -1 
NTOT=1C24 
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Nt;I^• = 23 
IFFIL=r 
K K = 10 
1 )=KK 
2) =': 

M(3)=r 

FACl=5.2r2/(.S^PHO^PSLOPE) 

FAC? = UMFAN/IJMVCLT 

FAC3=CABS(H(1)) 

rAC4 = FAC 1/ FAC3 

FACE=F AC2/FAC3 

FAC6 = FAC2/ (CA3 S( H( FRUN'+n ) ) 

C 

RAKGE( 1 ) =G . 'J5 
RANGF( 2}=< .3 
PANGE(3)=15.0 
RANGF(4)=-15.3 
C 

rr 1C^2 1=1,512 
FREQ(1)=I-1 

n02 TIMEITd )=(I-1 )«l.C/l''24. 

01 1 I=1,MBSKIP 

1 CALL IF'DATAJ DATA) 

L = 0 
C 

nC 11 SIDE=1,2 
DC 11 HC'LE = 1,25 
L = L+ 1 
C 

DC 3 IBLK=1,4 
J= ( I PLK-1 ) <'1024 
CALL IMDATA(OATA) 

DC 3 JJ=l, 1024 
1=J+JJ 

3 U( I )=CATA(JJ) 

C 

CALL HREAN(U,NCBLKS,BMU(SIDE,HOLE),BMSU(SICF,HOLE)) 
C 

CALL BZERC (iJ, K IN,KOUT ) 

C 

KPHASEI SIDE,mLE)=KOUT 
KSTART=KOUT 
KLAST=KSTART4-1023 
C 

PMJ(SICE,hCLE)=FAC5*RMU(SIDE, HOLE) 

FF SU(SIDE,HOLE) = (FAC6**2)*8MS0(SIDE,HOLE) 

C 

DC 6 K=KSTART,KLAST 
I=K-KSTAPT+1 
6 A( I, 1, 1)=FAC2<'U(K) 

C 

If SFT = -1 

CALL hAFR(A,M,INV,S,!FSET,IFERR) 

C 

A( 513, 1,1)=A(513,1,1)/H(512 ) 

A( 1, 1 , 1)=AI 1,1,1)/HI 1 ) 

Gxn) = 2.^*'A(l,l,l) ♦CUNJGI A( I, 1, 1 ) ) 

PC 21 1=2,512 
KCONJG=NTOT- I 2 
A( I,1,1)=A(I,1,1)/H( I) 

CCN=CCNJG( A( I , 1 , 1 ) ) 

A(KC(NJG,1,1)=CCN 

21 GX(I) = 2.0<‘A(l,l,l) «CON 
C 

SUM=C .C 
DC 22 1=2,511 

22 SLR=SUR+GX(I) 

50^=50^+. S*GX( 1)4-.5<'GX(512) 

C 

IFSFT= 1 

CALL HARM!A,M, INV,S,IFSET, I FERR ) 

DC 7 1=1,1024 
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c 


c 


7 ij( n =A (1,1,1) 


Vm R I T E ( 6,12 2 ) 

WRITE( 6, 1C 2) . ^ 

Vf P IT E ( 6, 1( A ) NRUN , S I DE f HOLE 
VsKITEI 6,12(. ) BI^U( SIDE,HOLE) 
CALI UTPLDTITIMEIT ,U,50,RANGE 


1 , 0 ) 


11 Ct'NTIKUE 


C 


C 


61 


5:' 1 

c 

i:-ei 


c 


WP ITE(6,1C C ) 

WPITF(6,lv';2) ^ ^ 

WPITE(6,lfA) NRUN, SIDE,HOLE 
WRITE(6,105) SUM 
Dl 61 LIKE = 1,128 
1X1= (L IME-1) ’i‘4+1 
f X2=IX 

lEIdlNE^EO.44).OR.(LINE.EQ.89)) WRITE(6,107) 
WPn E(t,K 6) ( IX,GX( IX), IX=IX1, 1X2) 


Wf- ITF ( 6, 119) 

WP ITE( 6, 10 2) 

WRITF(6,1C4) (NRUN,SIDE,HOLE) 
WPITE(6,12C) QMU,( S IDE ,HOLE ) 

DE RCCl LINE=1,89 
I X1=(LINE-1)«4 + l 


IX2=IX1+? 

IF(L.EC.44) WRITE(6,1C7) 
WRITF(6,121)(IX,U(IX),IX=IX1,IX2) 


rX' ICC 1 J=1,512 

GX( J)=GX(J)/SUM . 

CALL DRAW'(48' ,FRE0,GX,0,^,LArtC(L) 
C,r ,0,8,5,G,LAST) 

CALL PRAW(3),TIMEIT,U,0,2,LARE1, 
*•: ,C',f ,C , 8,4,0, LAST ) 


, IT1L4,6C .0 ,0.2 
ITIL2,r .004,0.0 


RANGE!3)=9.0 
RANGE!4)=-.R0 


, 


, 


C 

C 

c 

c 


c 

c 


L = 0 


DC 13 SIDE=1,2 
nC 1? Hf]lE = l,25 


L=L+1 


J= ! 1 RLK-1 ) >:^1C24 
CALL INrATA!DATA) 

DC 2 JJ= 1, 1124 
I =J+JJ 

P!I)=DATA!JJ) 

CALL EiMEAN ! P, NOBLKS , BMP ! S lOE, HOLE ) , 3MSP! SIDE, HOLE ) ) 

CALL EZEKC!P,KIN,KCUT) 

KSTAPT = KrUT«-l 
KLAST=(<STARTfl023 


C 

C 


C 


FAC7 = FAC1/!8MU!S I DE , HOLE )’<‘*2 ) 
HMP!SIDE,HCLE)=FAC4’«'BMP!SIDE, 


HOLE)/!HMU!SICE,HOLE )**2) 


VA 4 K=KSTART,KLAST 
1=K-KSTAPT+1 
4 A! I, 1, 1 ) = FAC7*P!K) 


CALL HARM! A,M, INV,S,IESET,IFERR ) 
A!513,1,1)=A!513,1,1)/(G!512)’'H!512)) 
A ( 1, 1, 1 ) =A ! 1 , 1 , 1 ) / ( G ! 1 ) ’»‘H ! 1 ) ) 

GX ! 1 )= 2 * A! 1, 1, 1 ) *CnNJG! A ! 1, 1, 1 ) ) 
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c 

PG IQ 1=2,512 

A ( I , 1, 1 )=A ( I , I , 1 ) / |G( I )''H( [ ) ) 

KCCNJC=NTCT-I+2 
CL N = CC’NJG( A( I , 1, 1 ) ) 

A(KCCKJG,1,1)=CaN 

19 GX(I ) = 2.r »A(I,1 , 1 )ACOG 
C 

SUK = C 

nr 2C 1=2,511 

20 SLM=SUN' + GX( I ) 

SGM= .5^GX ( 1 ) + . 5*GX (t)12 ) + SI)M 
C 

IF StfT= 1 

CALL HAPMI A,M, INV, S, IFSET, IFEP.P ) 

Ot 5 t=l,lC24 
P( I ) = A(I,1 , 1) 

5 P( I ) = P ( I ) 4-e,MP( S IOF,HOLE ) 

C 

PL 12 1 = 1,1^24 
IF=( L-1 )>iaC24 + I 

12 ^»PITE( 1* IB) P( I) 

C 

WF11r(6,119) 

UPITL(6,1C 1) 

WF ITFI 6, 1i.!4) riRUN, SIDE , HOLE 
UF ITE (6,12'') 3MP( S IOC,HOLE) 
nc 8rpr line = 1,39 
IX 1=(L INF-1)*4 + 1 
IX2=IXl+3 

IF(LINE.EC.44) URITE(6,107) 

I"'-'' wriTE(6,114)(IX,?( IX),IX=IX1,IX2) 

C 

DC per 2 1 = 1, 1024 
h:C 2 P( I )=P( I )-ei'^P( SIOF,HJLE) 

UFITE(6,122) 

UPITEI6,U 1) 

UPITF(6,1C4) NRUM,SIDE,HOLE 
WPITE(6,12C) BPP(S IDE,hole) 

CALL UTPLC'KTIMEIT ,P , B'', P ANGE , 1, ?) 

C 

Ur ITE(G,10C) 

WPITE(6,ir1) 

UP1TE(6,1l4) NRUN, sice,hole 
U hlTE(6,ir5) SUP 
DC 6f LINE=1,12P 
IX1=(LINE-1)*4+1 
I X 2= IX 1+3 

IF((LINF.FC.44).0H.(LINE.EQ.89)) WRITE(6,1C7) 

60 UFITE(6,1C6)(IX,GX(IX),IX=IX1,IX2) 

C 

PC’ P: 5 0 1=1,512 
r.5C GX( I ) = GX( I )/SUM 

CALL DPAV.(48C,FREQ.GX,C,'‘,LABC(L) , ITIL3,6C-.3,0.2, 
,C,C,C,8,4,>^,LAST) 

C 

13 CONTINUE 
C 

UKITE(2,5rC0) RWP 
C 

DC 5L IT = 1,1024 
no 51 H0LE=1,5' 

I r = ( K LF-l )*1''24+ I T 
PfADdMB) Y(HCLF) 

SinF=l 

IFOLF=HOLE 

If (HOLE.GT.25) I HOLE = HCLE-25 
ir(HC)l E.GT .25) Sir)E = 2 
51 Y(HOLE)=Y(HOLE I-bMP(SIDE, IHOLE ) 

C 

WP ITF(2,5P(3 ) Y 
C 
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c 

c 


DC. bl H0LF = 1 f 25 , ^ » 

Y(l-'OLF) = Y(Ht)LE + 25)-Y( HOLE ) 
b2 X (HOLE )=Y( HOLE)*ETMHOLE) 

CALL QTFG(ETA,Y,Y,NOIKJ 
CK(IT)=Y(NCIM) 

CALL OTFGIFTA,X,X,NDIM) 

FATA! IT) = y (NnI^^) 

50 CCNTIKUt 

SUH1=C .C 
SLN2=r .0 
DU 27 SIDE=l,2 

sf F'l-SL!°i +EMU (S lOE, HOLE) 

Slji^2 = SUN'2+RMSU(SIDE»HCLF) 

71 CCMINUr 

DC 36 SIDF=1,’ 

36 PMSutsKiF^HQLE) = SQRT( ^.OI^SQ^^.TIBMSUISIOE.HOLE) )/ 
*bN-L( SinEtHCLE) 

BFUA=SUMl/5'!. 0 
PK SUA=SUF'? /50 .C 
Of- A = BMUA 

RAT K;=SQPT (2 .0 *L)MS a )/UMA 
QyEAN=.5X-'RhOy^UWA«*2 

Y( HOLE )=GN'P( 2 'iHnLE) -BNPI 1 t HOLE ) 

29 X (tlUL F )=Y( HOLE I^FTAI HOLE ) 

CALL OTFGIFT A,Y,v,NCI P) 

FNKEAN=Y(MDI M) 

CALI QTFGI FTA, X,X,Nni“) 

CPPFAN=X(NCIM) 

n CATMI )=(CPMEAN+nATA( [ ) ) / I FNME AN + C N11 ) ) 

N r RI K ^ = 1 

CALL BMEAN(DATA,N06LKS,CPM,CPMS) 

t P IT R 61 i r ! N R UIM , F R U N , R A T 10,0 M A , A L P H A 
WFITElcIllZ) FNMEAN 
D(; 32 LIMF = lt89 
I CM = ( L II E-l ) '^A + 1 

IF (lIne'^^IfO. 4A) k^RITE(6,10 7) 

32 vJ FITE (6,113)1 lCNfCMICN)f ICN-1CN1» ICN2) 

RMTElbllK ) NRUNfFPUNtRATIOtUMA,ALPHA 

RPITF(6,111) CPP 
DP 33 LINE=1.59 

IC Fl = ( LlNF-1 )'-^4 + l 

IC P2=ICP1^3 

iMn"LnSi,icp,o5{IiicO?iCP=icPi.icp?. 

WR IT£(6,218) NRUN 

37 UFITE(6?214) * HOLE,ETA(HOLE)t(BMSU(SIDE.HOLE),SIDE- 

wPITE(6,2'''P)NRU'‘'l,UPA,OMSA,RATin,OMEAN 

WPITF(6,?13) NRUN 

28 WPlft(6?21A)’hCLE.ETAlHOLE),(BMP(S10E,H0LE). 

* S I PE = 1 , 2 ) 

WPltE(6,215) NRUN,FNMEAN 

DC 34 1=1,1324 

34 A( I , 1, 1 ) = CN( I) 

CALrHAPMl A,M. INV, S.IFSET.IFERR) 


, 2 ) 
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DC 35 1=1,512 

Gy{I)=2.':<^A(I,l,l) ^COKJGl A (1,1,1)) 

35 SI K=SUM4-Gy ( I ) 

SC> = SUN-.‘^ <'GK( 1)-.5''GX(512) 

rtF ITT ( h, F C ) 

V^FITE(6,115) 

V^FITf(f;,llC) NKUN , F RUN, RAT 10, DM A .ALPHA 

WE ITE(6, 1( 5) SUM 

DC ?P LIC:E=1,128 

ICM = ( LlFF-1 )’^4 + l 

IfN2=ICNl+3 

IF((L1NE.FC.44).0R.(LINE.FQ.3'^)) WR1TF(6,107) 

3b WF ITE( C, 1< 6) ( ICN,GX{ ICN) , ICIM= ICNl , ICN2 ) 

PC ICC 4 H0LF = 1,2 3 
1304 Y(H0LF)=-BMP(1,H0LF) 

call DFAW(23,FTA,Y,l,3,LAbEl, it IL5,3.2,1.C , 
^-0,C,C,C, 8,4,0, LAST) 
n( ICC 5 HC LF = l ,2 3 
ir5 Y(HrLE)=-bMR(2 .HOLE) 

CALL [)RAW(23,ETA,Y,3,5,LArtEl,ITIL5,''.2,l.:, 

«C ,C,a, 4,0, LAST ) 

CALL rRAU,(3'',TIMFIT,CN,^,2,LABEl,ITIL7,C.rC4,:. 

,C,C,8,^,0,LAST ) 

DC ICC 6 1=1,512 
i:C6 GV(I ) = GX( I )/SUM 

CALL PPAW ( 480 , FftFO.GX ,0,: , LABEL, IT IL6,6.'.0 ,. 2, 
,C,C,C,a,5, j, LA ST) 

DC ?P 1 = 1 , 1 .■'24 
39 A( I, 1, l)=|iATA( I) 

IESFT=-2 

CALL H APM( A, M, INV , S, I FSET , 1‘^ERR ) 

SUM=r.O 

DC 4f. 1=1,512 

GX ( I ) = 2.LX-A( I , 1 ,1 ) <=CONJG( A( I , 1 , 1 ) ) 

4 0 SL.M=SUM + GX ( I ) 

SLM=SUM-.5*GX(1)-.5*GX{512) 


WE ITP( 6,1C'' ) 

WE ITF(6,116) 

i^EITE(6,110) NRUN.FRUN,RATIO, UM A,alpha 

WFITE(6,1C5) SUM 

(1L 41 LINr = l ,128 

ICP1=(LIME-1)*4+l 

ICP2=ICPl+3 

ir((l INL.eC.44).nR.(LINF.EQ.8'P)) WRITF(6,107) 

4 1 WPITE(6,1C 6)( ICP,GX( irp) ,ICP=ICPI , ICP2) 

C 

DC ILCR 1=1,512 
n'8 GX{ I ) = GX( I )/SUM 

CALL CRAW(483,FREO.GX.O,',LABEL,ITIL8,6C.0,0.2, 
-XC ,C ,0,0,8,5,0, LAST ) 

C 

DC 42 1=1,1)24 

4 2 P ( I ) =CM I ) 

K K R = 9 

CALL RhAPM(P,KKR, INV,S, IFERR ) 

WP I TE ( 6, ICf' ) 

WPITt( 6, 117)(P( I ) ,I = l,n?6) 
r'( 4 3 1=1,10 24 

4 3 P( I )=DATA( I) 

CALI PHARM ( P ,KKR, I NV, S, IF'^RR ) 

WE IT F(6,ICC ) 

WR ITE ( 6,11 7) ( P ( 1 ) , I-=l , 1026) 

CALL PhAEM (iJ.KKR, INV, S, IFERR ) 

WPITF(6,i:C) 

WRITE!6,117)(U( I) , 1 = 1,1)26) 


, 
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f 


KO F( PMAT (in t///////fT5 3, 'TABLt' ,/♦ 

*TA6, 'POWER SPECTRAL DENSITY') ^^ i i 

m FORM AT(T42*'UNSTEADY PRESSURE COEFFICIENT ) 

FORMAT (TAG ,'UNSTEADY VELOCITY') . r a r,. r^. 

format(T 42 ♦'RUN' , I2»T53f'SlOE'fI 2fT66» TAP ,13) 
FC RM AT (TA8 ,'ME AN ‘^GyARE= ' , F7 . } ) 

FCKMAT ( 13X ,A( 'GX( ' .13 ,' = » ^ • ^ *tr hmt i mhp n i . /l 
FC'RM AT (IFl ,///////. TA7, ' TABLc (CONTINUED) ./) 

FORMAT(IHl 'TABLc* ,/. 

=!=TA 1, 'UNSTEADY NORMAL FORCE COEFF I C I ENT S ) 

1C9 FORMAT(IHl,///////.TE3t'TABLt',/. 

’^TX3,'UNSTEADY CENTER CF PRESSURt ) ^ 


1'2 
1 DA 
I 3 5 
1C 6 
1C 7 
1S8 


110 FCRMAT(T12,'RUN' ,I2»2X»'FREQUENCY=*»I2 * * ’ 

a ' E PS = ' . F5 . ,' ME AN VELOC ITY= ' . F6.2. ' FT/ScC ♦ 

^2X,'ALPHAS* ,12,' DEG' ,2X,'SF=n2A H2') 
format (TA2,'MEAN CENTFR OF 


111 
11 2 
113 
llA 

115 

116 
117 
1 IS 
120 
121 
122 
2'^ 8 


PRESSURE=',F6.3) 


FC PMAT (TB'^ ,' MEAN 
FORM AT(11X,A( 'CN( 

FC RMAT (11X,A( 'CP( ' , 13, ■ )=• ,n.^, 

FORMAT(TAl ,'UNSTEADY 
FORMAT(TA3,'UNSTEADY 


F('RMAT(1 

*T22,'W AS',F7.2,/, 
=)‘T32,'THE 
XT32, 'TEE R ATlU 


C0EFF1CIENT=',F6.3) 


NORMAL FORCE 
' ,13,' )=' ,F6.3,11X) ) 

,.r^,M = ',F5.3,llX)) ^ , . 

NCRMAL FORCE COEFFICIENT') 
CENTER OF PRESSURE') 

FORMAT ( IX , 1 )E1C . 2 ) ^ r-.. 

r-rRMATIU'l ,///////,T53,'TABLE' ) 

FCRMAT(T5^- ,'MEAN= ' ,F8.3) ,-,y i i 

FORMAT(1]X,4('U( ' , 13,') = ' ,F7.3,1DX ) ) 

FC R). AT ( IF i. ) ^ ^^ ^ /WErAGF MEAN VELOCITY FOP RUN',13, 
2 ,, , 

AVERAGE'mFAN square via S * , F7,2 , / , 

OF UNSTEADY TD STEADY 'WAS ' , F7, A, / , 
^T^2.'THE AVERAGE Q i^AS',F7.2) 

213 FORMAT!IHl,///////,T53,'TABLE' ,/, 

-T32,'MFAN pressure COEFFICIENTS',/, 


/, 


=« T 3 2 , ' R UN ' , I 3 , / , 


'CP UPPER' 


*^T32, 'TAP',T39,'LOCATION',T5' 

J(j(' 'CP LOW^P*) 


215 


RUN', 13 ,/ 


~ ■> . 


’^•T~2,'IS',F6.3) 

218 FLPMAT(1H1,///////,T53,^ABLE',/, 

^:-T32,'RATIO OF SOR T ( 2 ) =^5 CRT ( V ^ LOCI T Y 
«• MEAN SQUARE) TO MEAN VELOCITY',/, 
• Rl INi* , I V , / , 

f 


^T32,'RUN-, 


’!'T32,'TAP' ,T39,'L0CATICN' ,T5'‘ ,'UPPER' 


T6C,'LOWER') 


8'' 1 EL'RM AT ( 5 I 5 ) 

503 FCRMAT(5F1C.C) 
5': 2 FORM'AT ( 1 6F 5. S ) 
5:A FORMAT(8E1C.3) 
535 fGPMAT(?CAA) 

8 IDO FORMAT(BOA X) 

STCP 

END 


SLePrUTINE 0TFG(X,Y,Z,N0IM) 

DIMENSION X(1) ,Y( 1 ) ,Z(1) 

SUM2=C . 

DL 2 I=2,NDIM 

SUM2-SUM?+0.5MX( I ) - X ( I -1 ) )-x ( Y ( I )'^ Y ( I - 1) ) 
Z(1-1)=SUM1 
Z ( NDIM) = SUM2 
return 

END 
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sir ROUT I KF HMF AN( P, \nHLKS , 8M, ) 

h I It N P( 1 > 

MK=C 

3NS=: 

81^='' .C 

1 DC b I=1,NC3LKS 
J=(I-1 )^-l-24 

suM=c.e 

nr 3 11 =^ 1 , 1 '24 
JJ=J+II 

IF(^1K .FC. 1) GO TU 2 
SCt^^SUM + Pt JJ ) >**? 

Gr TC 3 

2 SUt''=SOM + P( JJ > 

3 CfNTIKUf 
SUP=SCN'/ir24 

IF (PK ,EQ. ') GO TO 4 

8FS=5PS+SUy 

Gf TO 5 

4 8P = BV+SUN' 

5 CrKTir.UF 

IF (MK .tC. ?) GC TO b 

B.''S=PPS/NCBLKS 

G( TC 8 

6 3M = FFV(M(.F L KS 
CO 7 I=1,JJ 

7 P( I )=P(I»-BM 
MK =1 

GC TC 1 
3 CCNTINOF 

I'l FfRVAT ( •MOnLKS= • , I 2,/, 2.:. X, • MEAN=' , FI 7.4,/,2C X, 

* = ' ,F1( .4,//) 

1:3 FCRKAT(2'X,'the FCLLCWING ARRAY LISTS THE INPUT TIMF H 
*,/, 2CX,'ZERr MEAN AND NORMALIZED .^ITH RESPECT TC TH 
’<-F • ,//) 

1C2 FfRKATI IX, nEll.3» 

PETtiPN 

END 


SLPROUTINF INDATAICATA) 
niRFKSILN OATAIKPA) 

1 FORMAT(16(64A4)) 
RCen(^,n DATA 
rftorn 
END 


St PFCUTINF RZEkfK X,KIN,KCIIT ) 

I I MF NS ION X( 1 > 

X 7 FF 0= I. r 
K0UT=1 

m 1 K = 21 ,KIN 
IF (X (K ) .GT . GO TO 1 

IF (X (K + n .LT. (• .0 ) GO TO 1 

I F ( ( X( K-n .GT. O.'M .OP. ( X (K-2 » . GT.C .I ) GO TO I 
CK=APS(X(K)I 
CKP1=APS(X(K+1I) 

IF(CK.LT.CKPII GO TO 2 

KT RV = K + 1 

XTRY=CKP1 

IF(XTRY.LT.X7FRO I GO TO 3 
G( TC 1 

2 KTPY=K 

X T R Y = C K 

IF (XTRY.LT .XZEPCI GO TO 3 
GC. TC 1 

3 KnUT=KTRY 
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F^ TRY LTHfL'K X,Y,NOATA,RANGE , KK Z , MOGCIJR ) 

444 IF (F-CnCUR. EQ.l.OR.MOOCUR.EQ.l > JSET = ') 

JSET=JSET+1 

IF(JSFT.GT.4) JSET = 1 

n( T '.r I =1 ,NOATA, KKZ 

I PTX = 6'(YMAX-Y(I))/YRANGE + 1.5 

IPTV=F.:.>;-( X( I )-XMIN)/XRANiGE+l .5 

IF ( IPTX.GT.61.OP.IPTY.GT.8l) GO TO TO 

IF ( I PTX. LF . ; .OP. I PTY. LE.C ) GO TO 7.> 

GP in(IPTX, lOTY ) = XCHAR(JSET) 

GL TC 7lG 
7C IFRR=IFRF+] 

7:o ccNTir-uF 
C 

C CC^^PUTF PROPER SCALE NUMBERS 

C 

IF <P.C['CUP. EO.T.rR.POnCUR. EO. I ) GO TO 3000 
I F (MCOCUR.EQ.2) RETURN 
GO TO 022 

aOE'') XlNCP=XRANGE/4. 

YIN(R=YRAMGE/6. 

XSCAlL(1)=XMAX 
YSCALF(I)=YMAX 
GO PC 1 = 2,5 

80 XSCALF(I)=XSCALE{I-l>-XINCR 
DO 81 1=2,7 

81 YSCALF(I ) = YSCALF( I-l)-YINCR 
C 

C OUTPOT SFCTICN ^.ITH GRAPH 

C 

IF(RCDCUP.EQ.O.CR.yCOCUR.EQ.3)G0 TC 022 
RETURN 

17 F( RMAT ( 12X , 

1 IP E1^.3,4(IC X,E1'^. 3)/15X,2H<'*,8(lPH + ^’l'**'i>***T 

722 RFITE(6,17) XSCALE(5) , XSC ALE(4),X SCALE(3),XSCALE(2),XS 

11=1 

I = " 

DC Kl IK=1,61 
IT ( I )«! ,ci,92 

91 VnFITE( 6,18) YSCALE ( I n , (GR I0( IK, IX ) , IX = 1,8 1 ) , YSCALEl I I 
13 FCRMAT(3X,1P 

1 Fl0.3,2X,lH+,lX,8lAl,lX,lH+,2X,Fn.3) 

II = I I 
GOTO K2 

92 RhITE(6,19) (GRID( IK , I X),I X = I , 81) 

19 FCRPAT(15X,IH*,IX,81A1,IX,IH* > 

102 1=1+1 

IF ( I-i:I 1C 1,103,1:3 
193 1=C 
101 CCNTINUF 

RRITF(6,22) XSCALE(5),XSCALE(4),XSCALE(3),XSCALE(2),XS 
22 FC RN AT ( 1 EX ,20*^^, 3( 10H+* = »** + *-** ) ,'<H + *'i'/ IP 

1 12X,EK.3,4( lOX 

IF(IEPR) ICOC, IC 1 

19' 1 WP ITE(6,2'' ) I ERR 

2' rCPPATd'CX 'NUPUER of points out of RANGE =' 14) 

i:‘0 rftupm 
fnc 
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OOO wO ooo ooo o oooo ooo 


ir>'V{ 1 )=',■ 

CL 98r l. = l,MT 
INV{ L'Me yp 4 i) = mtlexp 

DC 07r J=2,LM1FXP 
Jj=J+LPlEX F 

97 2 IfJV(JJ) = INV(J)+NTLEXP 
LEXf =KTLEXP/2 

9 8: IVIE XP = Lf^l EXP’^: 

7S2 If- ( U bET) 1 2,89 5,12 
E^ C 


SLHPCUTINf UTPLCT (X ,Y ,NDATA,RANGE ,KK2 , VQnCUR) 
C I MENS ICNi GR in { A1 , 81 ) , XSCALEf 5 ) , YSC ALE ( 7) 

01 YENS I CM X (1),Y (l),f<ANGr(4) 

IMEGE-K’t'2 G7I0,BLANK,r)CT,XCHAP{4)/lH. ,1H+, 1H-^,1HX/ 
DATA rUT , QL AMK/7.4B41 , 74"-^'' / 

GkID is the MATRIX USED TO PLOT THE POINTS 
IERP=" 

X*” AX = RAKGt ( 1 ) 

Xiv* IK = PAMGL ( 2 ) 

YPAX=RANGE(3) 

YK IN = h ANGE (4) 

CHECKITC X AND Y ^CINTS AMD PLOTTING THOS<= OUT CT 
AT THE PAPGIN 

Dfj 3.' 1 = 1 , EDAT A,KK,Z 

IE(X (I)-XMAX) ?.r,5,2:.5,22G 
■’2 A X ( I ) = XMAX 
IE RR = IERF +1 
Gt'TO 2r 

2‘5 IF {X ( I )-XM IN ) 2" 3, ? I', 2n 
2 ' 3 X ( I )=XMIN 
IFRR=IEFP+1 

31C If(Y (I)-YMax)215,215,?12 
21 2 Y (I )=YMAX 
I ERR= IEPP + 1 
GOTO ?: 

215 iriY ( I )-YMIN )217, 30,3'" 

217 Y ( I )=YMIM 
If HE =IEKP+ 1 

30 CE NT I HUE 

PLCiTTIMG X AND Y AXIS , IE NECESSARY 

XF AMCt =XFA X-XMIN 
YC AMGC = YVAX-YM IN 

BLANKING GUT MATHIX-IGRIO) 

DC 1 = 1,61 

DO ?C! JJ=1,B1 

ni GF ini I , JJ)=BLANK 
3'0 CC NTINUE 

Y TFST=YMAX*YWIN 
XT CST= XMAX*XMIN 
IF(XTEST) 1 ,222,222 

2 2 2 IE (YTE ST) 333,444,444 

1 lY AX I S = «'' . x-( - XP IN ) / XRANGE + 1 . S 
DC 4C 1=1,61 
4 0 GP ini I,IYAXIS)=OOT 
GCTC) 222 

3 3 3 I XAX I S=6. . X'YwAX/YKANGE+I . 5 

i)C 6C I=!,P1 
60 GPID( IXAXIS, I ) =00T 

PLACING POINTS IN THEIR PROOEP GRID PuSITIONs 


R ANb 
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GC TC 

921 IP1=INV( JPD/NTVNl 
931 I = 2^( 1PP1+ IP 1) + 1 

IP (J-I) 3^',945,94b 
940 T=A(I) 

A( n=A( J ) 

A(J)=T 
T=A( 1+ 1) 

A( I + 1)=A(J+1 ) 

A( J+n=T 
945 CGMINUr 
950 J=J+2 
86C JJ1=JJ1+JJC1 
C END CF JPPl A^JD JP2 

C 

870 JJ2=JJ2+JJC2 
C END CF JPP2 AND JP3 LCCPS 

C 

380 JJ3 = JJ3+JJ03 

C PNC CF JPP2 LOCP 

C 

990 IF( IPSET)841,B95,995 
391 U(. 89? I = 1 , NX 
992 A(2*1) = -A(2*1 ) 

9^5 RETURN 

TPE FOLLLNING PROGRAM COMPUTES THE SIN AND INV TABLES. 

iro M1=MAX1(K(1),M(2),P(3)) -2 
PT = PAXr (2,MT ) 

9C 4 IF (f^T-ID) 906 ,906,905 
90 5 IFEFP = 1 
GO TC 695 
9.06 IPEPR = C 
NT=2**PT 
NTV2=NT/2 

SET UP SIN TABLE 
THETA=PIE/2**(L+1) FOR L=1 
91C THFTA=.7853981634 

JSTEP=2**(PT-L+1) FOR L=1 
JSTEP=NT 

JDIF=2**(PT-L) FOR L=1 
JDIf=?**(MT-L) FOR L=1 
JD IF = NTV2 

S(JDIF)=SIN(THETA) 

DC 95F L=?,MT 
THETA=THETA/2. 

JSTEP2=JSTFP 
JSTEP=JOIF 
jr IF = J STEP/2 
S(JOIF)=SIN(THETA) 

JC l = NT-jniF 
S(JCl)=CCS(THETA) 

JLAST = NT-J STEP2 
IFIJLAST - JSTEP) 950,920,920 
920 DC 94C J=JSTEP,JLAST,JSTEP 
JC=NT-J 
JD=J+JDIF 

940 S(J0)=S(J)*S(JC1)*S(JDIF)*S(JC) 

950 CONTINUE 

SIT UP INV(J) TABLE 

960 PTLFXP=NTV2 
C 

C PTLFXP=2**(MT-L). FCR L=1 

LP irxp = i 
C 

C LP1FXP=2**(L-1). FCR L=1 
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C PIT-Pf VFFSFT. the FOLLOWING KGUTINE PUTS THEM IN OP'TF 

NT SU = I'.T’’=NT 
M2PT=K ^-NT 


c 

160 

I r (M3 

MT ) 

37 j, 36<* 

c 

36 3 

M3 GP 
1CC3 = 

. ( P 
1 

FJ. MT 


N'^ VNT=N?/NT 
MINN3=NT 
GL TO 3P2 
C 

C K3 LfSS TP!AN MT 

370 ICC'’=2 
N ? VNT=1 
i\TVN3=NT/N? 

MINN3=N3 

330 JJD3 = MTSC/N3 
M2MT=M2-MT 

•+50 IP I r'2MT )47 ; , 46' ,46' 
C 

C M2 GP. rp PK MT 

460 IGC2=1 

N2VNT=N2/NT 
M I N''i2=NT 
GC TO 480 
C 

C M2 LfSS THAN MT 

47C ICC2 = 2 
N2 Vf T= 1 
NTVr.2=NT/N2 
MINN2=N2 

430 Jjr2=NTSU/N3 
N'l f/T = M 1-MT 

S6C IF CM MT)670,550,560 
C 

C Ml GP. OP FO. MT 

560 IGC1=1 

N 1 VN T= N1 /r- T 
MINN1=NT 
GO TO 560 
C 

C Ml LESS THAN MT 

57C IGC1=2 
N 1VNT= 1 
NT VM = NT/M 
»''I NN1 = M 

5P0 JJD1 = NTSC/M 
OOO .JJ3=l 
J = 1 

no P5' JPP3=1,N3VNT 
IPP3=INV(JJ3) 

PC B7C JP3=1,MINN3 
Gf' TO ( 6 1r ,62'') , 1GC3 
Mv'' I P3= IMV( JP3 ) «N3VNT 
G( TO 6 3C 

62C IP 3=INV(JP 3)/NTVN3 
630 I3 = ( IPP'^+I P3)*N2 
7''0 JJ2=1 

DC 87 ’ JPP2=1,N2VNT 
IPP2=INV(JJ2)+I3 
0( 8 60 JP 2 = 1, MINN 2 

GO TO (7ir ,72:),IGC2 
MC I P2= INV( JP2I’!'N2VNT 
GO TO 730 

720 I P2=l MVUP2)/NTVN2 
7 30 1 2 = ( 1 PP2+I P2 )*M 
HC'' JJ1=1 

DO JPP!=l,MVNT 

1PP1=INV(JJl)+12 
Df' 8 5C JP1 = 1,MINN1 
GC TO (8U ,320 ) , IGCl 
310 IP 1= IMV( JP 1) ♦NIVNT 
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17C I?C=-I3C 

W3i 1)=-S(I3C) 

W? (2 ) = S(I3CC) 

GC Tn 

180 Vv3(l) = -1. 

Ui 3 ( 2 ) =. 

GC TC 2''C 

3’'^ I IN THIRO gUADRANT 
I3CCC=NT4-I3CC 
I3CC = -I3CC 
w3(1) = -?(I3CCC ) 

W3(2)=-S(I3CC) 

200 ILAST=IL+JJ 

nr 2?C I = JJ, ILAST, IDIF 

KlAST=KL+I 

C(; 22( K = I tKLAST , 2 

Kl=K+KBIT 

K2=K1+KF I T 

K3=K2 + Ke IT 

PC TWO STFPS WITH J NOT P 
A(K>=A(K) + A(K2)’:'W2 
A(K2)=A(K)-A(K2)^^2 
A( Kl ) = A( K 1 ) ’kW + AI K3 )*W3 
A(K3)= A(Kl )^W-A(K3 ) ^W3 

A(K)=A(K)+A(K1) 

A(K1 )=A(I' )-A(Kl ) 

A(K2) = A(K2)*-A(K3)=M 
A(K3) = A(K2 )-A( K3)-^I 

P= A(K2 »’^W2 (1 )-A( K2 + 1 ) ><W2 ( 2 ) 

T=A(K2)’^W2(2) + A(K2+1)>^^W2( 1) 

A(K2)=A(K)-R 
A(K)=A(K)+R 
A(K2+1)=A(K+1)-T 
A(K+1)=A(K4l )+T 

R=A(K3)*W3(l)-A(K3+l)*W3(2) 

T=A( K3 (2 ) + A(K3 + 1 ) ’I'W 3( 1 ) 

AWR=A(Kl)*W(l)-A(Kl4l)*W(2) 

AWI=A(Kl)=fW(2)+A(Kl + l)*W(l) 

A ( K? ) = AWR-R 
A(K34l)=AWI-T 
A ( K ! ) = AWR4R 
A(K141) = AWI4T 
T=A( Kl ) 

A( Kl) = A(K > -T 
A(K)=A(K)4T 
T=A(K 1 + 1 ) 

A(Kl+1 ) = A(K + 1 )-T 
A(K+ n =A(K + 1 ) + T 
R = -A(K3+1 ) 

T=A(K3 ) 

A( K3 > = A(K2I-R 
A(K?) = A(K2 )+R 
A( K3 + 1) = A(K2 + 1)-T 
220 A(K2+1)=A(K2+1)+T 

Ere OF I AND K LOOPS 

23G JJ=JjniF+JJ 

FNP OF J-LCOP 

235 Jl AST = A’F JLAST+3 
240 C( NT IMIF 

END OF L LOOP 

250 COM INUE 

Er 0 01 IP LOOP 

WE NOW HAVF THE COMPLEX FOURIER SUMS BUT THEIR ADDRESS 
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R (K?+l) 

T = A(K2) 

A(K2) = A(K)-R 
A(K) = A(K)4-K 
A(K2+1)=A(K+1)-T 
A(K^1 )=A(K4n + T 
C 

Atr,R = A( K1 )-A< Kl + 1 ) 

AW I = A<K1 + 1 )<-A(Kl ) 

R=-AJK3)-A<K3+1) 

T=A( K3)-A( K34-1 ) 

A(K?)=(AWR-R)/R00T2 
A(K3+1)=(AWI-r)/RCCT2 
A{ K 1 ) = ( AWR+R ) /ROOT2 
A(K1+1)=(AWI+T )/ROOT2 
T= A{K1 ) 

A(K1)=A(K)-T 
A(K)=A(K)+T 
T = A( KH-1 ) 

A(Kl + l ) = A( K+n-T 
A(K+1)=A(K4l ) + T 
R = -A(K3+ 1) 

T= A(K ? ) 

A{K3)=A(K2)-R 
A(K?)=A(K2 ) 4R 
A< K3+J )=A(K?+l )-T 
^5 A( K2+1) = A(K2 + 1 )+T 

ir(JLAST-l) 235,2i!5,«?: 

90 JJ= JJ + JJOIF 
C 

C NfW DC, THF REMAINING J*S 

Cl 23( J = 2,JLAST 

F r TC H W • S 

DF F- W = W*=i= INV( J ) , W2 = W<'*2t W3 = W=i‘*3 
96 I=INV(J41) 

98 IC=NT-I 

w ( u =s (ic;) 

W(2)=SII) 

12 = ?*'l 
I2C=NT-I2 

IF ( I 2C ) 12r ,l 10 ,l''0 
C 

C 2=*'I IS IN FIRST QUADRANT 

1 W2 < 1 ) = S( I2C) 

W2<2) = S( 12 ) 

GC TC 13r 
lie W2(! )=. . 

W2(2)=l. 

GC TP 13C 
C 

C ?=< I IS IN SECOND QUADRANT 

120 I2CC = 12C4NT 
I2C=-I2C 
W2(1)=-S(120 
W2(2)=S(12CC) 

130 13=1+12 

I3C=NT-I3 

IF (13C)16F,150,1^ 3 
C 

C 13 IN FIRST QUADRANT 

14C W3(1) = S(13C ) 

W^I2) = S<13 ) 

GC TC 20C 
15C W3(1)=C. 

W?(2)=l. 

GC TC 2rf. 

C 

160 1?CC=I3C+NT 

IF ( 13CC)19r,180,1TC 
C 

C 13 IM SFCCIND QUADRANT 
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Kt-' ir = KBIT/4 
KL=KFIT-? 

c 

c nc rnp J=‘: 

rc ar i = i, ili, idif 

KLAST=I+KL 
UC 3C K=I,KLAST»2 
K1=K+KBIT 
K2=K1+KBIT 
K ?=K 2+KfU T 

DC TWO STC.PS WITH J = r 
A ( K J =A (K ) +A( K2 1 
A(K2)=A(K)-A(K2I 
A( K1 ) = A( K1 H-A{ K3 ) 

A(K?) = A(K1 )-A(K3) 

A(K)=A(K)+A(K1) 

A( K1 ) = A(K ) -A( K1 ) 

A(K2 )=A(K2)+A{K3)*I 
A(K3)=A(K2)-A(K3)«I 

T=A(K?) 

A(K?)=A(K)-T 
A ( K ) =A(K ) + T 
T=A(K2+1) 

A( K2+ n=A( K+1 ) -T 
A(K + 1 )=A(K + 1)+T 

T=A(K3) 

A( K3 ) = A( K1 )-T 
A(K1)=A(K1)+T 
T = A( K3«-l ) 

A(K3 + 1) = A(KU1)-T 
A(Kl + 1 )=A(Kl + 1 ) +T 

T=A( K1 ) 

A(K1)=A(K)-T 
A ( K)=A(K) + T 
T=A(Kl + 1 ) 

A( K14-1 ) = A( K«-l ) -T 
A(K+1 )=A(K + 1 H-T 

R=-A(K3+1) 

T = A(K3) 

A(K3)= A(K2 )-P 
A(K2)=A(K2)+R 
A(K3 + 1 )=A(K2+1 )-T 
30 A (K2 + 1 ) = A( K2 + n +T 

ir (JLAST) 235,233,32 
32 JJ = JJP1F ♦■1 

DC FOR J=1 
ILAST= IL +JJ 
DC 35 1 = JJ, HAST, IDIF 

KLAST = KL+I 
00 35 K=I,KLAST,2 
K1 = K+KBIT 
K2 = Kl+KPIT 
K3 = K2 + K0 IT 

LF TTING W=(l+I )/RCCT2,W3=(-1 +1)/ROCT2,W2 = I , 
A(K)=A(K) + A(K2 )*I 
A(K?)=A(K) -A(K2KI 
A( K 1 ) = A( K 1 I^-W + AC K 3) ’’Wa 
A ( K3 )= A( K 1 )*W- At KU *W 3 

A(K)=A(K)*A(K1) 

At K1 ) = Atk)-AtK 1 ) 

AtK2» = AtK2 I + AtKa)-*-! 

AtK3) = AtK2)-At K3) -^I 
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OOOOD 


X2 ERC=XTPY 

IF (XZt Rt,.LT.O.CC:;''>l ) GO TO 5 
1 CLOT INUP 
5 RFTUPK 
Er.D 


SLOPrUTINE tURM( A, y, IPiV, S, IFSFT, IFERR) 

01 yFNSK.IM A( 1 ) , INV (1 ) , S( 1 ) ,N( 3 ) ,M( 3 > , NP( -^ ) ,W( 2 ) , W 2 ( 2 ) » 
EGLIVALEPXF (N 1 ,N( 1 )) ,{N 2 ,N( 2 )),(N 3 ,N( 3 )) 

10 IF( lAHS(IFSET) - 1 ) < 500 , 900,12 

12 yTT=yAxc (y (1), y( 2) ,N( 3)) -2 
P( 0X2 = SQRT( 2. » 

IF (yiT-MT ) lA,14,13 

13 irFFP=l 
RFTUPN 

14 IPEFR=^ 

Ml = y(1 ) 

P2 = M(2 ) 
y?=M(3) 

,\ji = 2*=«'yi 
\2 =?»X',V|2 
N3 = 2*’<‘M3 

16 IF(IFSET) 18 , 18 , 2 ' 

18 NX= M*Fi 2 >^N 3 
Fr = \y 

nr 1 ° I = i,Nx 

AIOX-I-I) = X(2*I-1)/FK 

19 A( 2 »I ) = -A( 2 * I ) / FN' 

20 NP ( 1 )=M*? 

NP(2 )= NP( 1) *N2 
NP(3)=NP(2)#N3 
C( 2 5r I D= 1, 3 
IL = ^P(3)-NP(ID) 

III = IL + 1 
PI = y(IP) 

IF (MI lOF" , 250 , 3 ' 

3 0 iriF=MP(lD) 

KP IT = NP(ID) 

M[ V = 2 ’'^ (M I / 2 ) 

IF (MI - MFV ) 6 C, 6 C, 4 '' 

C 

C y IS CDD. DO L =1 CASE 

40 KFIT=KBIT /2 
K(=KeiT -2 

DC 5 C 1 = 1 ,IL 1 ,IDIF 
KL AS 1 = KL 4 1 

DC 5 r K=I,KLAST ,2 
Kf =K 4 KDI T 

DC ONE STFP WITH 1 = 1 , 8 =' 

A(K)=A(K)+A(Kn) 

A( K(' ) = A( K ) -A(KD ) 

T=A(KD) 

A( KD) = A(K)-T 
A(K)=A(K)+T 
T=A(KD+l) 

A(KD +1 ) = A(K 4 l)-T 
50 A(K 4 1 )=A(K 41 ) 4 T 

IF (MI - 1 ) 250 , 250 , 5 ? 

52 LFIPST =3 

C OFF - JLAST = 2 >''*^(L- 2 ) -1 

JLAST= 1 
cr TC 7 C 

c 

C y IS EVEN 

60 LF IRST = 2 
JL AST =0 

70 DC ? 4 l L = L FIRST,M I ,2 
JJCIF=K 8 IT 
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TAPF CUNVERSUrj PROGRAM 


TH^ ASSERPLLR LARGUAGE PROGRAM 'FORM* THAT DOFS THE 
ACTUAL PIT CC'r vrRSinR WAS WRITTEN BY PIPCOPM C.ZELENY 
OF THE NAVAL POSTGRADUATE SCHOOL COMPUTER FACILITY STAFF 


C this PRC'GRA.m converts 7 TRACK, 24 BIT WORD TAPE 
C TO R TRACK 32 BIT WORD TAPE. 

DIPEN SION INDATAI1024) ,DATA(I'24) 

FACTOR = ir; . V( ) 

REWIND 2 
REWIMC 4 
NC BLK S = 1 

C DO J=1,K WHERE K IS THE NUMBER OF BLOCKS OF 1C24 EAC 
Dll 31 J=i,r:ro 
REAn(2,3,ERR=5C) INOAta 
CALL FCRM( INOATA) 

DO 1 1=1,1024 

1 DATA ( I )= INDATA ( I )-<rACTOR 
WRITE (4,^) DATA 
CALL P REAM (DATA, LMPLK S , 3M , PM S ) 

WR ITE(6,70 )J,BP,BMS 
GL TC 31 

5 ' WRITE (6,S 1 ) J 
31 CONTINUE 
: FCRPATI16(64A4)) 

51 F(RMAT (•^ ,5X,'READ ERROR, RECORD NC.=',I6) 

o6 F( RHATdX, I'EIC.2) 

71 F(U<MAT (ICX ,'record 00=',16,' ME A N= ' , Frt . 3 , • MEAN SO 

F 1 D. -^ ) 

STOP 
Ilf D 


SI f’KOUT I ME 3MFAN( p ,N0HLKS , BM, bMS ) 
DIVFE SION P( 1 ) 

ME=C 
BMS = ' .0 
BM = '' .f 

1 D(. 5 I=l,r;CBLKS 
J=( I-l 24 
SUM=0.6 

DO 3 I 1 = 1 , 1'24 
JJ=J + I I 

ir(MK .FO. D) GO TO 2 
SLV=SUM+E'( JJ )**2 
GO TO 3 

2 SUM=SUM+P(JJ) 

3 CONTINUE 
SUM=SUM/ K 2-^ 

I f ( MK . FD. :• ) GO TO 4 

BN S=FMS+SUN 
GO TO 5 

4 HN=PN+SUM 

5 CONTINUE 

IF ( MK .EC. "• ) GO TO 6 

PMS=BMS/NOPLKS 

Gi TG e 

6 PN=PM/NPPLKS 
Df 7 1=1,JJ 

7 P(I)=P(I)-RM 
XK =1 

GO TO 1 
6 CCNTINUL 
RETURN 
END 
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SURP;(niT ir F Pf)RM( INCATA ) 

FORM FTAPT ■' 

* SLIP FCUT INF FCP ^ ( I NOATA ) 

«' This SLf.FiniTINF 'aILL CCNVFRT 24 HIT RIN^RY WQRCS STORED IN 
X' -VN ARRAY LFNGTh SPECIFIFO dY THF INDEX VALUE TO 12 HIT PIN 

* \r,C PLACE TftCSE SAME WOPI')S HACK INTO INDATA 


* 

* 


LOOP 


CAT A 
NUM 


STM 14,12,12(13) This SUBROUTINE CONVERTS 


E/ LH 

6, ; 

24 

BIT 

BIf-ARY WCRUS TP 

USING 

X , f) 

12 

PI T 

w»)kri s 

U S I NG 
SE 

1 

DATA,7 

7, " 

11,=F »1C 24» 

THIS 

I S 

THE INDEX VALUE 


L 12,-' (1 ) 

L 2,NUM(12) 

LP 2,7 
SEPL 2,6 
SEL 2,2 
SE DL 2 , -> 

SPL 1,2 
SPPL 2,0 
SRL 2,2 
SPC'L 2,0 
ST 3,NUM(12) 

LA 12,4(12) 

BCT 11,LOOP 

L M 2,1 2, 2 » ( 1 3 ) 

MVI 1 2(13),X » FF* 

H(,P 15,14 

OSECT 

DS 1 F 

FT P 
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DISPLAY PROGRAM 


StBPCUTIKtS THAT GIVE THE POINT PLOT ON THE I3M 2250 
DISPLAY (JHIT V-tPE WRITTEN BY LT.(JG) G.J. VURHUFF, 
DEPARTMENT (F AF ROIM AUT I CAl FNGINEEklNGt 
NAVAL PLSTGRADUATE SCHOOL 


INTEGER HOL? 

DIME NS ION Y( G''-) tCPUIZb) ,CPL( 25 ) ,RAN'GE (4) , ETA (25 ) , 

^CEM(2,25) 

FFAL-^ TITLE!15)/'UNSTEACY PRESSURE COEFFICIENTS VS ET 

' / 

NDIM=23 
PANGL(1)=1. 

RANGE(2)=C . 

RANGE!? ) = 6 . 

RANGE!4)=-l. 

READ!5,1C 1)!ETA!I ), 1 = 1,25) 

1 1 ri'RMAT ! l6F5..r ) 

CALL COGSP 
R1R7 CC. NTINUE 

DC 1 IT=1, 1 24 
READ!^,1':2 ) Y 
K'2 PIRMAT!5CAA) 

IF!IT .NE.1 ) GO TO 3 
I =C 

DO 2 I=2,5C,2 
L = L+1 

CE L(L)=-Y! I-l) 

CPM! 1, L)=CPO(L ) 

CPL!L) =-Y! I) 

2 CPM! 2, L) =CPL< L ) 

GL TO 6 

3 nt A I =1,25 

A CPU! I ) =-Y! I I+CP*^! 1 , I ) 

DC 5 1=26,5 
J= 1-25 

5 CPL! J)=-Y! I )t-CPM! 2 ,J) 

6 C(NTINUF 

CALL LlTPfRT! ET A, C^U, MD IM , RANGE , 1, 1 , T I T LE , C R R 9 8 , CRR R<)) 
CALI UPH( R T! ET A , C L , NO I , R ANGE , 1,3 , T I TLE , IJ, I K , C9998 , C 
1 CC'NTIMUE 
9<^)8 CONTINUr 
CALL GGSP 
GC TO 9<-'97 
CONTINUF 
STOP 
EE C 


SUEPt UTI NE UTPCRT ! U , V , NDAT A , R ANGE ,K!<Z , MODCUR , T ITL F ,, * 

UTPCRT PrCCUCFS A POINT-PLOT CN THE 2255 GRAPHIC DISPLAY. 
U IS THE VECTOR OF ABSCISSAE 
V IS TtiL vrCTOC OF ORDINATES 

NDATA IS ThI length CF VECTORS U, AND V (THE NUMBER OF P 
range is TI'E SCALING VECTOR: 

RAEGEd) = MAXIMUM DESIRED X VALUE PLOTTED 

RANGE!2) = MINIMUM OESIREO X VALUE PLOTTED 

RANGE!?) = MAXIMUM OESIREO Y VALUE PLOTTED 

FAMGL!4) = MINIMUM DESIRED Y VALUE PLOTTED 

EVERY KKZ’TH ELEMENT OF THE INPUT LISTS WILL HE PLOTTED 
USER SETS MCDCUR = 

TITLF IS EITHER THE NA^E OF A 6?-BYrE ARRAY UR A 6;-CHARA 
CONST ANT. 

DATA SET FTC9F'''l MUST BE PROVIDED, nC3=! RECFM=F, LRE.CL=8C 
FIRST RETURN IS FOP RESTART, SFCONO FOR TERMINATION 


C 


CCMMLN/IGNORE/NULL , H^EV, IGS^NM 
niMENSIOF A!4),6!4) , X SC AL F ! 5 ) , Y SC AL E ! 5 ) 
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1 ,KULL ( 1 ) 19 ) , T ITLE( 15) 

OIMFKSIUN U( 1) ,V( 1) ,X( 1‘. ( ) ,/( I')'-) ,«ANGE(-i) 

PC AL>I 4 XCHAP < 4) / H+ , IH» , IHX , IH./ 

1 ^ Tf G1 k A, B 

I f ( MOnCUF . fo. • ) w'U Tfl(6,4 ) r I TLE 
4 FrPf'AT {•<■', 1 5A4/// ) 

CALL UTPLCiT (U , V , NO ATA , RANGE , KK Z , MOOCUR ) 

DC 2 I=l,f)DArA 
X( I )=U( I ) 

Y( I)=V ( I ) 

2 CfNTINUF 

IN'ITIALI/L GRAPHIC DATA SETS 


CALL I NODS(IDEV,IGCSNM) 

CALL SOATP(IGDSNM , 1 ) 

CALL SGDSL(IGDSNM,l.P, 1.64,1 I .3,1L .5 ,r .,0.,12.r I 2.) 
XMAX = RANbE ( 1 ) 

XMN=PANGr (2 ) 

YMAX=RANGF(3) 

VM IN=PANGE(4) 

XR AN(.L=X.VAX-XM IN 
YF AI GF = YMA X-VM I M 
OX =XPANCE/5S. 
nY=YP ANGE/^i, . 

XLL=XMIN-DX 
XIJR = XMAX+DX 
YLL=YMIN-PY 
VUR=YMAX+DY 

CALL SOATL(IGDSNM,XLL,YLL,XUR,YUR) 

CALL I N(.nS( IDEV, IGDSNN) 

CALL SLATM(IGDSNN ,3) 

CALL SDATl (IGDSNN , 0 ,, 7 3 , S 1 ) 

PLCT title 

CALL PTLXTdGDSNN , T I TL E , 6 ’ , NULL , NUL L , NU LL j'', 51 ) 
CALL EXF0( IGDSNN,C9998,£9990) 

PLCT GRAPH FRAME 


A( 1 ) = I 1 
B(1)=7 
A(2) = 1 1 
C(2)=^o 

A(3)=t9 
H(3)=49 
A ( 4)=o 9 
B(4)=7 

CALL STPnS(IGDSNN ,69,7) 

GALL PLINEIIGOSNN ,A,B,NULI,NULL,NULL,4) 

CALL FXEC(IGDSNN,£9998,£9999) 

CCMPlJTE PP( PEP SCALE NUMBERS 

8j-:0 X INCR=XPANGE/4. 

YINCR=YRANGE/4. 

XSCALE ( I ) = XyiAX 
YSCALE(1)=YMAX 
or ec 1=2,5 

3' XSCALF(I) = XSCALE( l-l)-XINCR 
DC B 1 1 = 2, 5 

81 YSCALE( I) = VSCALE( I-l)-YIMCR 

PLlTTING X AND Y AXIS, IF NECFSSARY 

YTEST=YMAX <YMIN 
XTEST=XMAX«XMIN 
IE (XTEST) 1 ,222,222 
222 IE(YTEST)333,444,444 
I CCNI I NUE 

CALI STPCSdGDSNM ,r.,YMIN) 
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CALL PLINfcdGOSNM ,"'.,YMAX) 

CALL EXFO( IGOSNM,C9Q98,C9Q‘5^) 

Gf.TO 2 22 
133 CCNTINUE 

CALL STPnSdGDSNM ,xyiN,0.) 

CALL PLINECIGDSNM ,XMAX,C.) 

CALL FXFQ( IGDSNM, &9994,C9999 ) 

444 CGNTIKUF 
C 

C PLLT SCALE hUMSERS 
C 

REWIND 9 

DC ir2 1=1,5 

WRITE (<^, IP )YSCALE:( I ) 

18 El’RHAKGl' . 3, • - • ) 

19? CCMIMJE 

wPITr(9,54) 

54 Ff PVAT(12X,•I • ,4( I3X, • I ' ) ) 

WRITE(a,22)XSCALE(5), XSCALE(4),XSCALE(3),XSCALE(2),XS 
22 Et.PE'AT (4y , 5( 4X,G1 ^.3) ) 

Ei\C FILE 9 
RfWIND 9 
C(’ IFl 1 = 1,3 

REAr(‘^,53)(AREA(K>,K=l,19) 

53 F(PPAT(1CA4) 

L INE = 4 3-( I -1 ) '■='1: 

CALL PTEXT ( IGOSN'N , AR E A, 12 , NULL , NUL L , NUL L , 5 , L I NE ) 

CmLI E XFC( IGDSNN,r.999e,S990Q) 

I'l CCMINUF 

DC: IPS 1 = 1,2 

RFAPIo,53) (AREA(K) ,K=1,19) 

LIKF=P-I 

CALL PTEXKIGOSNN , ARE A , 74 , NULL , NUL L , NUL L , 3 , L I N E ) 

CALL EXEO( IGDSNN,(19993,C9999) 

133 rCMlNUF 

ITFNPl=IGOSNM 
ITEP:P2 = IGDSNN 
If RR=C 
GC ID 5 
C 

C (NCTF: ENTRY UPHCRT ALLOWS FOR AQPING CURVES TG GRAPH WH 
C HAS Al READY BEEN LAID OUT -- CALLED BY PRNCRT) 

C 

entry UPHCPTd),V,NCATA, range,KKZ,MCDCUR,TITLE,IGOSNP, I 
C 

IF (MCDCLIR. E 1.3 )WR I TE( 6,4 ) title 

CALL UTHOLDI U, V,NI)ATA,RANGE ,KKZ,MnDCUP ) 

C(J 3 I=1,NCATA 
X( I )=U( I ) 

Y( I)=V(I) 

3 CFNTINUF 
5 CCNTINUE 

IGCSNP = ITE MPl 
ICDSr;N=ITEPP2 
C 

C CHECKING X AMD Y POINTS AND PLOTTING THOSE OUT OF RANGE A 
C 

Of ?l 1=1 , NDAT A,KKZ 
IE(X (I)-XMAX) 2J5,?-:5,220 
22'^ X ( I ) =XMAX +OX 
IFRE = I FRP+ 1 
G( TO 21( 

2'-5 IE(X ( I )-XNIN )2'3,210,2i: 

2^3 X (I)=XMIN -DX 
IFPP= IFER+ 1 

219 ir(Y ( I )-YMAX )213,2 15,212 
212 Y (I)=YMAX +DY 
IFRR= lEPP + l 
Gt TO 3' 

’15 If (Y ( I )-YNIN)217 , 30,3 3 

217 Y (I)=YMIN -OY 
IFRF = irRK+l 
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CCNTIMfl- 

C 

C PLACING PCINTS IN THEIR PROPER POSITIONS 
C 

If (PI OCUP. FO.^’.CK.PGnCUR.FO. 1 ) JSFT = '' 

JSET =JSFT + 1 
If (JSET.GT .4) JSET = 1 
0( 91 1=1 , N'lATA.KKZ 

RXCriAK=XCHAR(JSFT) 

xi=xn ) 

YI=Y( I ) 

CALL PTEXT(1GDSMM,RXCHAR,1,NULL,NULL,NULL,XI,YI) 
CALL F XCQ( IGDSNK,&999P,&9 999) 

PI CONTINUE 

Ir(PCPCUR.tJ.■.CO.POOCUR.FO.X)GO TO 922 
RE TUEN 

T2? CONTINUE 
C 

IF ( IFOR) 1C K* 'll'''' 1 
ICTl CFAiTINUE 
l-Lt^lNn c 

KP ITEI9,2t ) IPRR 

2' F(RPAT( 'NUMBER OF POINTS OUT OF RANGE = ',14) 
ENC FILE 
REWIND 9 

F LAI. (0,52,) (^PEA(K) ,K=1,19) 

CALL PTEXT ( 1 GiJSNN , AR E A, 4C , MULL , NHL L , NUL L , p , 4) 

CALL FXCU(IGDSNN,09998,G9999) 

DC'* C(.NTINUE 

C IF (PLTCUF.EU.3IKETUPN 

C 

C WAIT FUR FN(i-KEY ATTENTION 
C 

CALL t PEATN( I ATTN) 

CALL RQATT( I ATTN, I COOL,NULL,G99 94,£9999) 

CALL NP\TL(I ATTN,-1) 

CALL F F'ATL ( I ATTM) 

CALL T PGCS ( IGDSN^- ) 

CALL TPGDS(IGDSNN) 

RETURf 

9994 CENT IEUF 
RETURN 1 
9999 CLNTINUE 
RETUFN 2 
END 


SUePCUTIME PRNCRT ( X , Y , NN , MODC UR , T I TL E ,♦,="•) 

C 

C PRNCRT PRGDUCES A PCINT-PLGT ON THE 229'' DISPLAY, IT CAN 
C SEVFP4I CURVES ON THE SAME GRAPH, DEPENDING ON THE CALLED 
C mcDCUR: 

C ='■ , THIS IS THE ONLY CURVE FOR GRAPH 

C =1, TF-IS IS THE riRST JF SEVERAL CURVES FOR GRAPH 

C =2, THIS IS AN INTERMEDIATE CURVE FlJR GRAPH 

C =2, THIS IS THE LAST OF MANY CURVES FOR GRAPH 

C X IS THE VECTOR OF ABSCISSAE 
C Y IS TF^f. VECTOR OF ORDINATFS 

C NN IS THE LFNGTFI OF VECTORS X AND Y (THE NU*^RER OF PCINTS 
C TITLF IS FITHFR THF' NAMF OF A S''-HYTE ARRAY OR A ftP-CHARA 
C CONSTANT. 

C PRNCPT CREATES THE PROPER RANGE VECTOR AND CALLS UTPCRT 
C PCINTS OF TF(E FIRST CURVE ON A GF^APH. SUCCEEDING CURVES 
C ARC TRUNCATED TO THE SA^E RANGE AND PLOTTED WITH CIFFEREN 
C first return is FOR RESTART, SECOND FOR TERMINATION 
C DATA SFT FT^QF "'1 MUST EE PROVIDED, DC 8= ( RE CF M= F , LR F CL = B'' 
C 

Cf MMOr./I GNOME/NULL, IDEV, IGSPNM 
C 

0 1 Ml NS ION XtI), Y(51), RANGE(4) 

UrULLd ) ,APEA( 19) ,TirLE( 15) 

ECUIVALFME < R ANG^ ( 1 ) , XM AX ) , ( R ANGE ( 2 ) , X M I N ) , ( RANGE ( T ) 
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1 


(f^ANGE(4) ,YM1N) 


BRANCH CN ^TDCUR 

IF (MCDCUF . Eg."-.CR.,'^GOCUR. EQ. 1 ) GU TT 40 
Gl TC A'-' 

GAUGE IMFUT DATA AND CCYpuTE RANGE FOR FIRST OR ONLY CURV 
XF'AX=-l.F2r 


xrIN=1 .E2^ 
YMAX^- 1.E2C 
YNIN=1 .E2C 

lie 1 I=1,NN 


IF (X( I )-XMAX) 
XMAX=X(I) 

YX‘'AX = Y( I ) 

6 , , 2 

IF (X( I I-XMIN) 
XMIN=X(I) 
YXMIN=Y(I) 

3,3,7 

IF(Y(I)-YMAX) 
YMAX=Y(I) 

X\MAX=X(I) 

8,8,4 

IF(Y( I )-YMIN) 
YM IN = Y( I ) 

XYR I N=X( I ) 

CC NT IF OF 

5,5,1 


CALL LTHCRT PROPERLY 


4':'' ir (KfCCUR. E0.2.CR.MGOCUR.EQ.0) GO TO “S'-’ 

CALL UTPCR T( X , Y,NN ,RANGE i 1 , MUDCUK , T I T L E , 119998 , C9999 ) 

GL TO ?(" 

3: CALL UPHCRTIX,Y,FA ,RANGE,1,MCDCUR,TITLE,IGC,IG01,C9998 
70 I r-(PLC'CUR .NE . 3 IRFTURN 


LIST TRUrCATICA RANGE AFTER END OF MULTI-CURVE PLOT 


' 1 
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6''' 


REV, I NR <0 

RP ITE( 9,1' 1 ) YMAX , 
FI'PMATI '(^AX Y = • 

1 • N n Y = • 

UR ITE ( 9, ir<") XKAX, 
FCPMAT ( 'MAX X = ' 

1 'MIN X = • 

ENC FILE 9 
I CCSNM= IGf 1 
RCUIN'D 9 
DU 6C 1=40,52 
READI'^ ,<L I ) 

FORMAT(1GA4) 

J= I 


XYMAX,YMIN,XYMI N 
,G14.7,• AT X = • ,G14.7/ 
,G14.7, • AT X = SGIA.T) 
YXMAX, XMIN, YXMIN 
,G14.7,• AT Y = ',014.7/ 
,G14.7,• AT Y = *,014.7) 


(APEAIK),K=1,19) 


CALL PTLXTIIGOSNM, area, 74,NULL,NULL,1,0,J) 
CALL E XEQ( I0DSNM,&999S,C9999,f.62) 

CENTINUC 


UAIT FOR END-KEY ATTENTION 
CALL CREATNI lATTN ) 

CALL RQATTI I ATTN, I CODE , NULL , (19993 , C9999 ) 
CALL MPATL(IATTN,-1) 

CALI ENATL(IATTN) 

62 CCNTINUE 

CALL TMGOSUGO) 

CALL TMGDS(IGDl) 

RF TURN 
9< 58 CCNTINUE 
RETURN 1 
9909 CCNTIA.UF 
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P[TUFN 2 
END 


SbBPLUTI^F ()eLAY(N) 

CALSF TIFT LAPSE CF N HUNDREDTHS CF A SECOND (REAL TIME) 

IFTEGlP TIPF.DT 
TIMt = I TIKE(D ) 

1 CCNTINUE 

OT=ITIKt(0 >-TIme 
IFIOT.LT.NISO TU I 
Rf TUPF. 

FKP 


SI PRC'UTINE EXEQI IGCSN^I,*, 

COFTRCL DISPLAY REGENERATION — HOLD FURTHER GRAPHIC DATA 
BUFFER IS PFACv 

CLMKR'N// TTENT/IATNL 

CCKKON/IGNCRE/NULL,IDFV,IGSPNM 

CALL P CATNIIATNL , ICCDE , 1 ,NULL,:U ,32 ) 

IF(ICCDF.f C.-l)CALL START!&RTM8) 

I r ( I CODr . E C. 31 )CALL F I N IS H ( CAug 8 , GD'JQD ) 

If I ICl DE.FC.32)PETUPN 3 
CALL EXECIIGDSNM) 

CALL fFLAY(5) 

RETURN 

■ ClTTINUf 
RFTUfN 1 
5-)3 COF.TINUF 
PFTURN 2 
END 


SUBROUTINE GOGSP 

INITIAIIZF CSP, 2 25'*, AND SEVERAL WOR K-OA T A-S E F S 

CFKKCN/ATTtNT/IATNL 

CC KMflN/ IGM’RE / NULL , I DF V , I GSPN‘^ 

NULL=-5 

CALL INGSP ( IGSPNM,nijlL ) 

GC Tf'. 1 

cNTRY PFTVIDES FOR P E - I N I T I AL I Z AT I ON 

ENTRY GGSP 
CALL TKDFV(IDEV) 

I CfTTINUE 

CALL INDEV(IGSPNM,22,IDEV) 

CALL SALPK(ICEV) 

CALL KLITS(inEV,3) 

CALL r RATL(I DEV, I ATNL ) 

CALL ENATNd ArNL,'>,3l ,32) 

RETURN 

ERF 


SLOP OUT INF START 

RESPLNSF TO RE-IN I TI AL I ZAT ION COMMAND DURING EXECUTION 

CCKKPN/IGNORE/NULL, IDEV, IGSPNM 
C 

CAlL SALPM(IDEV) 
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CALL TKCFV(IOEV) 

CALL INOFV(IGSPNM,22fIDEV) 

CALL INITl^St IGCPGL,3,2,''.,0. , 1 . , 1 . , ^. , ?. , 4 8. , 34. , . TRUE 
DC 1 1=1,1C 

CALL PTEXTdGDPGL,•RESTARTING G SP . . . ' , 17 , NUL L , NUL L , NUL 
1 CCNTINUF 

CALL EXECIIGnPGL) 

CALL nELAYd'''?) 

RETURN 1 
END 


SLBRnUTIMF FINISH!*,’*') 

PESPCNSE TO TEKMINATICN PEGUEST DURING EXECUTION 
(ALLUV) CNF CHANCE TO CHANGE MINO) 

CCMMDN/IGNCKE/NtLL,IDEV,IGSPNM 

CIPFNSICN FNT(13) 

lata fvt/•(••1********** EXECUTION ABORTED 8Y USER *** 
CALL SALRN(IOEV) 

CALL TNDEV(IOEV) 

CALL INDFV(IGSPNM,22,IOEV) 

CALL I MTOS! I G DPGL , , 2 , . ,0 . , 1. , 1. , 3 . , 0 . , 48. , 34. , . TRUE 

INTENSIFY MESSAGE 

CALL DGSEQ(IGOPGLI 
DC 2 1=1,1C 

CALL PTFXT!IGOPGL,*LAST CHANCE... START OR FINISH?',31 
INULl ,F,17) 

2 CONTINUE 

CALL ENSEW(IGOPGL ,Krv ) 

CALL E XEC! IGUPGL) 

CALL CRATL(IDEV,IATNL) 

CALL FNATN!IATNL,',31) 

flash FNKB lights and CRT DISPLAY 

1 CCNTINUE 

CALL DELAY (5^) 

CALL CHIT! IGDPGL,NULL,KFY) 

CALL MLITS!IATNL,2) 

CALL DEL AY(3C ) 

CALL INCL!IGOPGL,NULL,KEY) 

CALL ML ITS(IATNL,3) 

CALL PCATNd ATNL , IC0DE,1,NULL ,C,-35) 

IF! ICODF.EC'. })GC TO 1 
CALL ML ITS ! I ATNL , 1 ) 

IF!ICCOE.EC.-l )CALL START(&Q9Q8) 

CALL RESET! IGOPGL ) 

CCNFIRM TERMINATION 

DC 3 1=1,1C 

CALL PTEXT!IGDPGL,FMT!2),44,NULL,NULL,NULL, 3,17) 

3 CLNTINUF 
CALL EXEC! IGOPGL) 

WF I TF(6 ,FMT) 

CALL OELAYIS'^O) 

CALL TMOEVIIOEV) 

RFTUPN 2 

9P93 CONTINUE 
PETUPN 1 
END 
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SLBRrUTIKF IN 1 TCS( IGDSNM,IOATM,ICH,X I,YI,X2,Y?,X3♦Y3 
INTFFFACC TC GSP GHS-CRF AT lOM A.MO QPT 1 JN-OEF I NI T I ON ROUTI 
CCPMGr./IGNORE/NLLL , lOF Vf IGSPNM 
LCGICAL LP 

CALL INGDS(IDEVflGDSNM) 

CALL SGDSL(I GOSNM,X1,Y1,X2,Y2,0., , 1. , 1. ) 

CALL SOATL(IG0SNM,X3,Y3,X4,Y4) 

CALL SDATMCIGDSNM,IDATM) 

CALL SCHAM(IGOSNM, ICHAM) 

IF(LP)CALL SLPAT(IGOSNM,n 
RF-TUHr’ 

Eh C 


SLOROUTir.E CREATN( 1 ATTN) 

INTFFFACE TC GSP 'CRATL* 

MAINTAIN RESTART/ARORT LEVEL AS ACTIVE 

rr P^CN/IGNCRE/NULL , IDEV, IGSPN'^ 

CALI fRATL(lUEV, I ATTN ) 

CALL RLPFC(I ATTN, I ,2,1 ) 

CALL NLPEOII ATTN,2,4,1 ) 

CALL PPATLdATTN, 1) 

CALL ENATNnATTN,G,31 ,32) 

RETURN 

END 


StfiROUTINE RQATTI I AT TN, I CODE , / I RR AY/ , «• ,* ) 
INTERFACE TC GSP 'POATN* 

RESTART OR ABORT ATTENTIONS TAKE PRECEDENCE 

CCNKCN/AT TENT/1ATNL 
COMMON/IGNORE/MULL,IDEV,IGSPNM 
C 

DIMENSION IRRAY(l) 

CALL RCATNI I ATNL, I CODE, 1,NULL ,'^,31,32) 

Ir(ICODE.EC.-l )CALL START(E99P3) 

1 F( ICODF.FQ.31 )CALL F1NISH(E999B,&9999) 

IF ( IC(J0F.E0.32 )RFTURN 
CALI MPATL ( I ATTN,-1 ) 

CALL ROATNII ATTN,ICODF,2 , I RRAY , ?,-35) 

CALL MPATL(1ATNL,-1) 
in ICOOE.CC.-l )CALL START(CQ998) 

IF( ICCOE.FO.3 1 )CALL FINISH(E999R,09999) 

RE TURN 

9993 CCNTINUF 
RETURN 1 
9999 CONTINUE 
RETURN 2 
FNC 
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